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ABSTRACT 

NGC 205, a close satellite of the M31 galaxy, is our nearest example of a dwarf elliptical galaxy. 
Photometric and kinematic observations strongly suggest that NGC 205 is currently undergoing tidal 
distortion as a result of its interaction with M31. Despite earlier attempts, the orbit and progenitor 
properties of NGC 205 are not well known. In this paper, we present an optimized search for these 
unknowns by combining a genetic algorithm with restricted iV-body simulations of the interaction. 
This approach, coupled with photometric and kinematic observations as constraints, allows for an 
effective exploration of a 10-dimensional parameter space. We represent the gravitational potential of 
M31 as a static analytic bulge-disk-halo model. NGC 205 is modeled as a static Hernquist potential 
with embedded mass-less test particles that serve as tracers of surface brightness. We explore three 
distinct, initially stable configurations of test particles: a cold rotating disk (no velocity dispersion), 
a warm rotating disk (small amount of velocity dispersion), and a hot, pressure-supported spheroid 
(isotropic velocity distribution with no rotation). Each of these models is able to reproduce some, but 
not all, of the observed features of NGC 205. This leads us to speculate that a rotating progenitor with 
substantial pressure support could match all of the observables. Furthermore, plausible combinations 
of mass and scale length for the pressure-supported spheroid model of the progenitor of NGC 205 
reproduce the observed velocity dispersion profile. For all three models, we find that NGC 205's line- 
of-sight distance and proper motion are well constrained by NGC 205 's kinematic profile and surface 
brightness distribution. Orbits that best match the observables place the satellite 11 ± 9 kpc behind 
M31 moving at very large velocities: 300 — 500 km s^^ on primarily radial orbits. Given that the 
observed radial component is only 54 km s~^, this implies a large tangential motion for NGC 205, 
moving from the north-west towards the south-east, that translates into a predicted proper motion of 
~ 0.1 mas yr~^. These results suggest that NGC 205 is not associated with the stellar arc observed to 
the north of M31 and to the north-east of NGC 205. Furthermore, NGC 205's velocity appears to be 
near or greater than its escape velocity, signifying that the satellite is likely on its first M31 passage. 
Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — galaxies: individual 
(NGC 205) — galaxies: interactions — methods: A^-body simulations — meth- 
ods: numerical 



1. INTRODUCTION 

In the hierarchical model of galaxy formation, low 
mass dwarf-sized galaxies collapse and merge to form 
the massive galactic structures that we observe to- 
day. Minor mergers in the present-day Universe 
continue to influenc e the properties of stellar halos 
(jBullock fc JohnstonI [2005), and perhaps the appcar- 
ance of galactic disks (IGilmore etal.l 120021 : lAbadi et aP 
l2003t lYoachim fc DalcantonI l2005l ). How this process 
of minor merging works in detail is not well under- 
stood. The Local Group is a prime locale for study- 
ing minor merging events. In the Milky Way, three 
satellites are undergoing tidal disruption and accretion: 
the Magellanic Stream (|Mathewson. Clearv. fc Murravl 



(M31), our nearest large neighbor, also hosts a 
number of tidally distorted and disrupted satellites. 
These include the dwa rf ellip tical galaxy NGC 205 
JKent 19 871: IChoi. Guh athak urta. fc JohnstonI |2002| : 
[Geha et all 2006f). the compac t elliptical galaxy M32 
(iKentI Il987t iChoi et al.l [200l. the Giant Southern 
Stream (llbata et al.l 120011 : iFont et al.ir2006l: iFardal et all 



Il974t IConnors. Kawata fc GibsonI |2006[). Sagittarius 
dwar f dibata. Gilmore. fc IrwinI 11994 iMaiewski et al.l 
120031 iNewberg et al.l I2003D. and Canis Major dwarf 



(jMartin et al.ir2005f ). The Andromeda Spiral Galaxy 
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|2007|). and possibly th e satellite galaxy Andromeda VIII 
( Morrison et al.l [20031 ). In this paper, we present a de- 
tailed study of the tidal interaction between M31 and the 
dwarf elliptical galaxy NGC 205. 

In addition to being building blocks of larger galaxies, 
dwarf galaxies are interesting in their own right. Of the 
various dwarf galaxy types, dwarf elliptical (dE) galax- 
ies are the least understood, and their origin is heavily 
debated. What is known about dEs is that they are a 
galaxy population found exclusively in denser regions, ac- 
counting f or more than 75% of the ob j ects in cluster envi- 
ronm ents (jBinnev fc Tremaind Il987t iTrentham fc Tullvl 
|2002| ). They are typically dark matter poor in their in- 
ner regions and predominantly gas poor. The two gen- 
eral dE formation theories postulate that these galaxies 
are primordial building blocks or an evolved/transformed 
population, but neither theory is capable of explaning 
all the observables. The first theory, based on Cold 



Howley et al. 



Dark Matter models, asserts that these are low mass, 
primordial building blocks that formed early (z > 6) 
in their present locale. If dEs formed by this sce- 
nario, we could gain insight into the very first galax- 
ies to form by studying these relics. The second the- 
ory postulates that dEs are born as more massive spi- 
ral galaxies which are subsequently altered as a result 
of interactions with other galaxies through tidal eff ects 
(IMoore. Lake, fc Katd [19981 IRicottifc Gnedin|[2005h . If 
dEs are transformed spirals then this transformation is 
expected to cause an incr ease in the velocity dispersion 
and induce some rotation (jMoore et al.lll99a l. However, 
this scenario is unable to explain the observed distribu- 
tion of rotation speeds and velocity dispersions among 
Virgo Cluster dEs. Specifically, it cannot explain the 
presence of a substantial population of dEs with no de- 
tectable rotation and relatively low velocity dispe rsion 
(|Geha. Guhathakurta. fc van der Mareill2002ll2003D . To 
complicate matters further, the faint underlyi ng disks are 
found in both rotating a nd non-rotating dEs ("Gch a et al.l 
I2003t iLisker. Grebel. fc Binggeh et al. 2006). It is also 
possible that the dEs result from some combination of 
the above two mechanisms. A complete model explain- 
ing dE formation is still needed. Local group dE galaxies 
offer an excellent resource for the exploration of these and 
other formation scenarios. 

Unlike the Milky Way which does not host any 
dE satellites, M31 hosts a total of three: NGC 185, 
NGC 147, and NGC 205. NGC 185 and NGC 147 lie 
far enough away from M31 to escape present tidal dis- 
tortion. NGC 205 lies a projected 37' from M31 and is 
our nearest example of a tidally distorted dE galaxy. It 
is a gas poor, low luminosity, early-ty pe galaxy with a n 
exponential surface brightness profile (jChoi et ani2002f ). 
defining features of dE galaxies. The little gas and 
dust in NGC 205 is concentrated within a 1' radius, be- 
yond which the satellite is essentially gas and du st-free 
^Young fc LolfT997HV^^Ich. Saee. fc Mitchel]|ll998D . The 
current line-of-sight distance estimate between NGC 205 
and M31 is about 39 kpc, wit h M31 residing 785 ±25 kp c 
and NGC 205 at 824 ± 27 (jMcConnachie et al.l[2005h . 
Throughout this paper, we adopt this distance of 824 
kpc to NGC 205. The absolute F-band magnitude of 
NGC 205 is My = —16. 5, corresponding to a few per - 
cent of M31's lumin osity ('Hodgd ll992HChoi et al.|[200l ). 
iBender. Paauet fcl ^ieto (1991) reports a mass-to-light 
ratio for NGC 205 of (M/L) « 7. This value is con- 
sistent with dynamical models of dEs that suggest a 
global mass-to-light ratio of 3 < Fy < 6 and is indica- 
tive of an intermediate to old stellar population contain- 
ing little to no d ar k matter within an effe ctive radius 
(jGehaet al.ll2002D . iDe Riicke etHI (|2006f) finds a B- 
band mass-to-light ratio within 2i?e = 260" (1.04 kpc) 
of {M/L)b == i.5tl]l{M/L)Q. This supports the con- 
clusion that there is little dark matter contained within 
the inner regions of NGC 205. It is not known whether 
NGC 205 resides in an extended dark matter halo. 

There is strong e vidence tha t NG C 205 is tidally 
distorted by M31. IChoi et all (|2002f ) measure twist- 
ing of elliptical isophotes, along with a subtle break 
in the surface brightness profile, providing confirma- 
tion of tidal distortion found in earlier studies. A ra- 
dial velocity survey of red giant branch stars (RGB) 



in NGC 205 by iGeha et al.l (|2006t ) notes an abrupt 
turnover in the semi-major axis velocity pr ofile coinci- 
dent with the tidal distorti on rep orted by IChoi et al] 
(|2002f) . iMcConnachie et al.l (|2004f) postulates that the 
~ 1° (15 kpc) long, arc- like feature seen in the northwest 
quadrant of M 31 could be a stellar stream a.ssocia ted 
with NGC 205. iDemers. Battinelh. fc Letartd (|2003f) re- 
port possible evidence of tidal debris to the west of 
NGC 205 beyond its tidal radius in their study of C stars 
in the system. Within an effective radius, NGC 205 ap- 
pears to be undistorted. However, beyond this radius, 
there is strong evidence for distortion. 

In this paper we investigate the large parameter space 
defining the orbit and history of NGC 205 using a re- 
stricted A^-body code and a genetic algorithm for opti- 
mization. Modeling satellite/parent galaxy interactions 
is extremely difficult when there are no gaseous or stellar 
streams to directly constrain the path of the satellite's 
orbit. In the absence of tidal trails, as is the case for NGC 
205, there is a limited set of observable quantities that 
can be used as constraints: stellar radial velocities, pro- 
jected sky position, and constraints on the line-of-sight 
distance. To complicate matters further, large param- 
eter spaces with inherent degeneracies typically define 
galaxy encounters. Hence, in weakly distorted systems, 
it is difficult to find unique, qualified fits to the observ- 
ables. For this reason, few optimized A^-body studies of 
mildly distorted galaxies have been carried out. In our 
simulations we employ a genetic algorithm to achieve op- 
timization (jHollandlll975l : lGoldberg|ll989( ). This particu- 
lar optimization technique has been successfully used in 
similar studies of mo r e heavily distorted dwarf galaxie s 
(jTheis fc Kohlel[20?)ll : iTheis. Gerds. fc Spinnekeii [200lh . 
Using a restricted A^-body code, we explore three dis- 
tinct, initially stable configurations of 1000 test particles 
representing NGC 205: a (cold) rotating disk without 
velocity dispersion, a (warm) rotating disk with velocity 
dispersion, and a (hot) spheroid supported completely by 
randomly distributed, isotropic velocities. Detailed pho- 
tometric and kinematic observations enable this unique 
study b y providing stron g constraint s on the simulated 
system '(|Choi et al.ll2002t iGeha et al.ll2006D . 

The paper is laid out as follows: in §[2] we discuss the 
photometric and kinematic observational data used to 
constrain our system. In §[3] we give the details of the 
parameters, galaxy models, and methods used in our nu- 
merical simulations. In §|4]we present the results of the 
numerical simulations. In §[5] we discuss the implications 
of our results. Finally, in §[6]we discuss the observational 
and computational improvements that could be applied 
to this study. 

2. OBSERVATIONAL DATA 

Our study of NGC 205 is unique because we enforce 
both photometric and kinematic constraints on our sim- 
ulated system. Although this cannot ensure that we con- 
verge on the correct orbital history of NGC 205, it does 
remove a significant amount of the degeneracy inherent 
in studies that only use photometric constraints. Strong 
evidence of tidal distortion is seen both in the photome- 
try and the kinematics of NGC 205 at a major-axis dis- 
tance of r ^ 4.5' - 5' (1.1 - 1.2 kpc) (jChoi et al.ll27)0l 
iGeha et al.ll2006f). A dis cussion of the surface photome- 
trv bv IChoi et al.l (|2002f) is included in §0 and a dis- 
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cussion of the internal kinematics by iGeha at al.1 ()2006f ) 
is included in S I2.2I These observations are summarized 
in Figure [T] 

2.1. Photometric Data 

The surface photometry of NGC 2 5 in t he B- and 
/-bands was obtained by IChoi et aD ()2002f ) using the 
Kitt Peak National Observatory 0.9/0.6 meter Burrell 
Schmidt telesc ope. Due to the in completeness of the I- 
band coverage, IChoi et all ()2002f ) employed the 5-band 
data and IRAF's ELLIPSE task to determine the sur- 
face brightness contours of NGC 205 out to a limiting 
magnitude of ^b = 27 mag arcsec"^, coinciding with a 
s emi-major a xis distance of 720"(2.88 kpc). 

IChoi et all (J2002) reports that NGC 205 is weU fit by 
two different exponential profiles at various semi-major 
axis radii. One profile, covering a semi-major axis range 
of 75" < r < 250" (0.30 < r < 1.00 kpc), has a disk 
scale length of r^^^ — 150" (0.60 kpc), where tq is de- 
fined as the radius at which the intensity has decreased 
by a factor of e; the other profile, covering a semi-major 
axis range of 150" < r < 250" (0.60 < r < 1.00 kpc) 
and r > 500" (2.00 kpc), has a disk scale length of 
Vq^^ = 170" (0.68 kpc). Also reported was a "subtle 
downward break" at r = 300" (1.20 kpc), which corre- 
lates nicely with the semi-major axis distance of r = 260" 
(1.04 kpc) that marks the location where the isophotal 
position angle and ellipticity stop increasing and begin 
decreasing, leading to a S'-shaped semi-major axis pro- 
file for NGC 205. The location of this isophotal twist- 
ing matches the radial velocit y turnover r a dius of 4.5' 
(270" or 1.08 kpc) observed bv lGehaet all (|2006D . thus 
supporting the conclusion that NGC 205 is interacting 
tidally with M31. 

We perform our own exponential fit of the form: 

I{r) = /oe-'-/'-" (1) 

to the intermediate parts of IChoi et al] (120021 ) inten- 
sity profile in order to estimate the disk scale length of 
the galaxy if it were tidally undistorted. We exclude 
the nucleus (interior to 100", or 0.40 kpc) where the 
data follows neither an exponential law (exp[— r/rp]) nor 
a de Vaucouleurs law (exp[— A:r^''"']), and the outer re- 
gions (exterior to 672", or 2.68 kpc) where the surface 
brightness and sky become comparable. We measure a 
disk scale length of rg = 148" (0.59 kpc) and a central 
(nucleus-free) surface brightness of /q — 21.1 magnitudes 
per square arcsecond. These results are used as a basis 
for generating an undistorted NGC 205 and modeling 
the inner and outer regions of NGC 205 in its final, in- 
tegrated state. Hence, we are able to remove the nu- 
cleus from NGC 205 by extrapolating the exponential 
surface brightness profile from the intermediate region 
(100" - 672") to the cent ral region (0" - 1 00"). The 
surface brightness profile of IChoi et al.l ()2002D . our expo- 
nential surface brightness fit, and the data used in our 
simulations are summarized in Figure [2] 

2.2. Kinematic Data 

A r ecent kinematic study of NGC 205 by iGeha et al.l 
(|2006) using the Kcck/DEIMOS multisht spectrograph 
resulted in the collection of radial velocities for 725 RGB 
stars. The stellar spectra were obtained using 4 masks, 



each covering an area of « 16' x 4' (3.83 x 0.96 kpc). Two 
of these masks were centered on NGC 205 along the ma- 
jor axis, while the other two were placed off center along 
the tidally distorted major axis, one in south-east (SE) 
and one in north-west (NW). These observations extend 
out to 20' (4.79 kpc), well beyond the tidal radius of 
NGC 205, and are spatially well distributed. The con- 
tamination by M31 field stars is estimated to be only a 
few percent. 

The systemic v elocity of NGC 205, as measured by 
IGeha etal.l (|2006D . is -246 ± 1km s~\ while the sys- 
temic velocity of M31 is — 300ib4km s~^ (|de Vaucouleurs! 
\imi ). This results in a velocity for NGC 205 of 
54 ± 5 km s^i relative to M31. The velocities of the NW 
tail of the satellite are more positive, while the SE tail, 
which is closer to M31, are more negative as compared 
to NGC 205 's systemic velocity. The opposite is true in 
t he central re g ions ( < 270"). 

iGeha et all (|2006( ) maps the semi-major axis velocity 
profile of NGC 205 by binning the RGB stars perpen- 
dicular to the S'-shaped major axis into > 1' (0.24 kpc) 
radial bins. The velocity of each bin, containing a min- 
imum of 25 stars, is determined by fitting a Gaussian 
profile to its stars. The resulting fit is plotted in Figure 
|3l The major-axis rotation curve of NGC 205 shows a 
distinct velocity turnover (with i^max = 9 ± 4 km s""'^) 
at 4.5', which is coincident with the estimated tidal ra- 
dius (~ 4') and o nset of isophotal twisting (see ? l2.ip . 
iGeha et all (|2006[ ) suggests that this feature is due to 
gravitational interaction between NGC 205 and M31, 
and infers that NGC 205 is on a prog rade encounter. 
If correct, over half of iGeha et al.l (|2006f ) evenly sampled 
RGB stars lie outside NGC 205's tidal radius, meaning 
they have been tidally stripped from the satellite and are 
no longer bound. The ratio of the maximum rotational 
velocity to the average velocity dispersion of NGC 205 
(assu ming an average e llipticity of e = 0.43) is Umax/c = 
0.21 (|Geha et al.ll200a) . Given that the expected ratio 
for an isotropic, oblate, rotational-flattened galaxy with 
e = 0.43 is vo/ao « 1.1 and that similar anisotropic 
Virgo dE's on average have wq/o'o < 0.1, NGC 205 falls 
somewhere betwee n a rotational supported body and an 
aniso tropic object (jBinnev fc Tremaiii3ll987l : lGeha et al.l 
l2006h . 

We perform a linear fit to the the interior regions 
(< 270") of NGC 205 whe re significari t tidal distortion 
has not yet occurred. Using iGeha et al] (|2006( l error bars 
for each velocity measurement, we determine a maximum 
and minimum slope for the satellite's central velocity pro- 
file of 0.0 km s~^ arcsec"^ and —0.08 km s~^ arcsec"^, 
respectively. This significance level corresponds to 1.83cr 
and is plotted in Figure |31 We use both this subset and 
the full dataset to constrain the motion of test particles 
in our numerical simulations (see S I3.6.2I and S I3.6.3P . 

3. NUMERICAL SIMULATIONS 

The previous section discusses the photometric and 
kinematic observations of NGC 205. In order to re- 
produce the observed features of NGC 205, we trace 
NGC 205 's center of mass back half an orbit, set up an 
initial configuration of test particles (to represent the 
undistorted NGC 205), and allow it to interact under 
the infiuence of both NGC 205's and M31's potential. 
The final, integrated galaxy is then compared to pho- 
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Fig. 1. — Observed properties of the M31-NGC 205 sy stem. The colo red dots represent NGC 205's semi-major axis velocity relative to 
its systemic velocity of —246 km s~^ as measured by Ge ha et al.l II2006I ). Note the turnove r in ve locity at the semi-major axis distance of 
~ 270" (1.08 kpc) is coincident with the isophotal twisti ng radius m easured bv lChoi et al.l II2002I ). The labels A-F are defined as follows: 
A=center of NGC 205, B=radius at which the use Choi et al.l 11200 21) data begins (internal to this we use an expon ential surface brig htness 
profile), C=radius at which isophotal twisting and radial velocity turnover occurs, D=radius at which the use of IChoi et al.l ||2002D data 
ends (external to this we use an exponential law), E=ccnter of M31, F=center of M32. T he projected distance betw een NGC 205 and 
M31 (points A and E) is 36.6', or 8.8 kpc if we take NGC 205 to be at a distance of 824 kpc (IMcConnachie et al.l|2005l ). The photometric 
contours B, C, & D lie at semi-major axis distances of 100" (0.40 kpc), 260" (1.04 kpc), and 672" (2.68 kpc), respectively. (Image courtesy 
of Phil Choi) 



tometric and kinematic observations. Strict upper and 
lower limits are placed on the variable parameters and 
a genetic algorithm is employed to effectively search our 
parameter space. In the following sections we discuss the 
explorable parameters of the system (iJ IS.ip . the poten- 
tials of both galaxies (S I3.2|) . the initial configurations of 
NGC 205 's test particles fS I3.3p . the integration scheme 
(§12131), the genetic algorithm used to improve on the ini- 
tial guesses (§ 13. Sp . and, finally, our method of comparing 
the integrated system to observations (ii l3.6p . The goal of 
these simulations is to determine the orbital parameters 
of NGC 205. 

For our simulations, the units are chosen such that the 
distance unit is 2503" (or 10 kpc if we assume NGC 205 
hes at a distance of 824 kpc), the velocity unit is 280 
km s~^, and G = 1. This roughly results in a mass unit 
of 1.8 X 10" Mq and a time unit of 35 Myr. 



3.1. System Parameters and Constraints 

There are a total of 10 free parameters in our numeri- 
cal simulation. Six of these parameters define the initial 
conditions of the undistorted NGC 205 system, while the 
other four define the final conditions of the present day 
observed system. The six initial condition parameters 
are the mass of NGC 205 (Af205), the disk scale length 
(ro.i), the relative orientation (in (p, 9) of NGC 205's ini- 
tial, undistorted disk, whether NGC 205 is on a prograde 
or retrograde orbit about M31 (A6'), and the Hernquist 
(1990) scale length (the radius at which a quarter mass 
is enclosed) of NGC 205's potential (a205). The four final 
constraints are the three-component final velocity vector 
of NGC 205's center of mass relative to M31's center of 
mass {vx,Vy,Vz), and the line-of-sight distance between 
NGC 205 and M31 (z). A detailed discussion of the 
upper and lower limits placed on each of these param- 
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Fig. 2. — Surface brightness profile of NGC 205 in magnit udes per square a rcsecond versus semi-major axis distance. The solid line 
illustrates the brightness profile i n B magnitudes m easured bv lChoi et al.l 1)20021 '). the dashed line represents the exponential fit performed 
using the intermediate regions of lChoi et al.l (|2002l ) dat a (100" to 672"), and the blue line defines the data used in our simulations. Note 
that in the intermediate regions (100" to 672") we use IChoi et al.l | |2002! ) measurements, and in the inner (> 100") and outer (< 672") 
regions we use our exponential fit. This effectively removes the nucleus from the inner region and models the data in the outer region where 
the surface brightness of the galaxy and the sky brightness become comparable. 



eters is motivated and discussed in detail below. The 
parameters are discretized into 256 steps within their 
given constraints, with the sole exception of the pro- 
grade/retrograde parameter, A0, which can take on a 
value of either or 1. This number of steps is chosen such 
that the genetic algorithm can thoroughly explore our de- 
fined parameter space, while limiting the time needed to 
perform the simulations. These given parameter values 
define approximately 10^^ possible orbits. In the contin- 
uous real world, there are an infinite number of orbits. 
A summary of the constraints is given in Table [T] 

Mass. The mass of NGC 205 is not currently well 
known. A recent constraint was placed on the mass con- 
tained within a 2i?e radius sphere by ' De Riicke et al.l 
(Ho p) of M205(2-R e ) = lO-St^.n x 10® Mg. We use 
De iDe Riicke et al.l (|2006[ ) minimum of Afaos > 8.2 x 
lO^Mgas a lower limit on the total mass of NGC 205 
(i.e. if all the mass were contained within 2i?e) and set 
a very liberal upper limit of A/205 < 5.0 x 10^ Mq. 

Disk scale length. This parameter is used when 
NGC 205 is modeled as a disk. In § 12.11 we determined 
that the observable disk scale length of NGC 205 is 
vq — 148". We assume that at earlier times NGC 205's 
surface brightness profile followed an exponential law 
governed by a different tq, which has since been affected 
by its interaction with M31. We therefore allow the ge- 
netic algorithm to explore the range 15" < ro,,; < 313" for 



the initial disk scale length in order to model the undis- 
torted NGC 205. If NGC 205 indeed hes at a distance 
of 824 kpc, then this tq range corresponds to 0.0625 kpc 
< ^0,1 < 1-25 kpc. The disk scale length is also used 
to determine the maximum distance an initial test par- 
ticle can be placed in the disk of NGC 205. We chose a 
maximum distance of 8ro,i; however, at this distance we 
expect to find very few particles since the initial distri- 
bution is subject to follow an exponential law. 

Initial disk orientation. These orientation parame- 
ters are used when NGC 205 is modeled as a disk. We 
assume axisymmetry for NGC 205 's initial, undistorted 
disk and describe its orientation in 3-D space by two 
of the three Euler angles. We define the plane of the 
sky as the x-y plane centered on M31 with the positive 
X-axis in the direction of increasing RA, the positive z- 
axis pointing towards the observed NGC 205 and the 
negative z-axis pointing towards us. In our defined co- 
ordinate system, the Euler angle i/i describes turning the 
galaxy's a:-axis clockwise in the x-y plane and is also al- 
lowed to vary from 0° to 180°. The inclination angle 6 
describes tilting the galaxy about the rotated x-axis and 
is allowed to vary from 0° to 180°. An inclination angle 
of 90° means the initial galaxy appears edge on, while 
an inclination angle of either 0° or 180° means the initial 
galaxy appears face on. Note the visual symmetry about 
180°, which does not translate into dynamical symmetry. 



Howley et al. 






iDD 



_ Dl 



I -- 



I 
1 

V 

^ -JD 



D — 



}1 



•1 



f 



i}****^f 



i*iJ 



i 



I 



il * } 



IW 




■1 


jlli*^ 


Kpj_-iiii_i- 


TO 


■^-■■■_. 



E -in - 
■S :n 

R' -^ - 



-lixf -MM} D im an TCC 

■ l-rilJMlilil_ 



4- 



-I- 



i 



iJ.. 



{ 



* jM H 



-LDtn 



S Bai - BrI ^'r An ■ ^E^m rd ]^ 



llQ] 



Fig. 3. — NGC 205's semi-major axis radial velocity (relative to M31) and dispersion profile as reported bv lGeha et al.1 1)20061 ) . (a) The 
black dots represent the combined velocity measurements of RGB stars. At distances of ~ 270" (1.08 kpc) from the center of NGC 205 
tidal distortion becomes evident, (b) The slope interior to the distortion of NGC 205. The ranges used to constrain our simulations are 
illustrated in the insert and correspond to a lower limit of —0.08 km s^^ arcsec"^ and an upper limit slope of 0.0 km s~^ arcsec"^ (or —20 
km s~^ kpc~^ and km s~^ kpc^^, respectively, when NGC 205 lies at a distance of 824 kpc). (c) The velocity dispersion of NGC 205 
along the semi-major axis. 



TABLE 1 
NGC 205 Genes 



Parameter 


Symbol 


Lower Limit 


Upper Limit 


Mass 


M205 


8.2 X 10« Mq 


5.0 X 10« Mq 


Disk Scale Length (disk only): 








...angular size 


ro,i 


15" 


313" 


...at distance of 824 kpc 


ro,i 


0.06 kpc 


1.25 kpc 


Clockwise Rotation (disk only) 





0° 


180° 


Inclination (disk only) 


e 


-180° 


180° 


Prograde vs. Retrograde (disk only) 


e-Ae 


A6» = 0° 


A6» = 180° 


Hernquist Scale Length: 








...angular size 


C205 


0" 


2503" 


...at distance of 824 kpc 


0206 


kpc 


10 kpc 


Distance from M31 


Z 


2 kpc 


76 kpc 


Center of Mass x- Velocity 


Vx 


-500 km s-l 


500 km s-l 


Center of Mass y- Velocity 


Vy 


-500 km s"l 


500 km s~i 


Center of Mass z- Velocity 


Vz 


49 km s-i 


59 km s-l 



Note. — Constrained parameters (genes) searched by the genetic algorithm in 
the numerical simulations. Each parameter contains 256 steps between the upper 
and lower limits, with the exception of 9 — A9 which can take on a value of or 
1. 
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The Euler angle ip describes revolution about the rota- 
tion axis and is meaningless due to the axisymmetry of 
NGC 205's disk. Similarly, the axisymmetry of NGC 205 
makes visual searches in the range 180° to 360° for (p and 
9 redundant. Thus, the position (but not the dynamics) 
of NGC 205 's initial disk is completely described by the 
parameters p and 9 in the range 0° or 180°. 

Prograde vs. Retrograde. This parameter is used 
when NGC 205 is modeled as a disk. Since the initial 
sense of rotation of NGC 205 is unknown, it is impor- 
tant to explore both prograde and retrograde orbits. A 
prograde orbit means that the galaxy is rotating in the 
direction of the encounter, while a retrograde orbit means 
it is rotating in the direction against the encounter. As 
mentioned in the previous paragraph, exploring 9 in the 
range of 0° to 180° provides a complete visual search, but 
not a complete dynamical search. For completeness, we 
introduce the rotation parameter A0 into our simulation 
whose purpose is to determine whether the Euler angle 9 
will be flipped about 180°. The rotation parameter can 
be set to either 0° or 180°, where A0 — 0° means 9 is 
unchanged and M = 180° changes 9 to 6*- 180°. Hence, 
this parameter essentially allows 9 to vary from —180° to 
180° so that both prograde and retrograde orbits can be 
explored. The motivation for using this separate parame- 
ter, A0, rather than just allowing 9 to vary initially from 
— 180° to 180° is that A0's binary nature provides an easy 
means for the genetic algorithm to instantly change the 
direction of the encounter via mutation or reproduction 
(see ? l3.5l for more information on the genetic algorithm). 
That is, the genetic algorithm is capable of changing the 
rotation of the disk by merely turning this parameter off 
or on. Note, this effect could have also been achieved 
using the 4> parameter instead. 

Hernquist scale length. We model NGC 205's po- 
tential as a spherical symmetric distribution of matter 
governed by a Hernquist profile fii l3.2[ Eqn[2]). The con- 
centration of the combined dark and baryonic matter in 
NGC 205 is governed by Hernquist scale length parame- 
ter, 0205 (or flo in EqnU). We set a lower limit on 0205 of 
0", corresponding to a Keplerian potential, and an up- 
per limit of 2503", which at a distance of 824 kpc would 
correspond to a scale length of 10 kpc. 

Velocity. The mean radial velocities of M31 and 
NGC 205 are -300 ± 4 km s"^ Cdc Va ucouleiui [l99lb 
and —246 ± 1 km s~^ ()Geha et al. 2006*), respectively, 
resulting in a relative radial velocity for NGC 205 of 54 
±5 km s~^ with respect to M31 (see %^^. Thus the 
radial velocity component, u^, is given a lower limit of 49 
km s~^ and an upper limit of 59 km s~^. The transverse 
velocity components, v^ and Vy, are given upper and 
lower limits dictated by the escape velocity of NGC 205's 
center of mass relative to M31, and are therefore not only 
functions of z and Wz, but also of M31's mass. In order 
to explore both bound and unbound orbits, we allow Vx 
and Vy to explore the full range of NGC 205's escape ve- 
locity, with a maximum escape velocity of 500 km s~^ if 
NGC 205's line-of-sight distance and velocity {z and v^) 
are at their minimums. 

Line-of-sight distance. The reported distances to 
M31 and NG C 205 are 785 ± 25 kpc a nd 824 ± 27 kpc, 
respectively (jMcConnachie et al.l 120041 ) . By adding the 
errors in quadrature, we approximate the line-of-sight 



distance to be 39±37 kpc and allow the genetic algorithm 
to search within the range of 2 kpc to 76 kpc for the 
final position of NGC 205. We recognize that this is a 
generous range to search, but decide to err on the side 
conservatism. 

3.2. Galaxy Potential Models 

The close proximity of M31 to NGC 205 stipulates that 
fairly accurate potentials are necessary in modeling the 
interaction between the galaxies. However, there is an 
implied limit to the level of complexity of the models re- 
sulting from the large number of orbits searched by the 
genetic algorithm. That is, using too simple of a model 
can result in incorrect final integrated positions and ve- 
locities, and, if the model is too complex, the genetic 
algorithm will spend large amounts of time attempting 
to complete a single orbit. Hence it is essential to find 
balance between the two requirements in order to come 
up with an optimal model. In this subsection we dis- 
cuss the potential models used in our simulations. We 
assume that the galaxy potentials are not changing sig- 
nificantly during the course of the interaction and thus 
remain static. 

Previous studies suggest dwarf elliptical galaxies do not 
contain a significant amount of dark matter in their in- 
ner regions. The amount of dark matter in the outer 
regions of dEs is completely unknown. Thus, one of the 
goals of this project is to probe the dark matter content of 
NGC 205. We make the simple, but appropriate, approx- 
imation that the potential of NGC 205 is a spherically 
symmetric spheroid with a Hernquist profile (Hernquist 
1990). However, it should be noted that the ideal model 
for NGC 205 is a stellar population orbiting in a fully 
consistent potential with a brightness profile lying some- 
where in between an exponential disk and exponential 
bulge (i.e. a hot, exponential disk). The difficulty of 
constructing such a model is demonstrated by the com- 
plex photometry discussed in § 12.11 For the purpose of 
this paper we model NGC 205's potential as, 

GM 
*6M = — (2) 

OQ + r 

where M is the total mass of NGC 205 (stars, interstellar 
gas, dark matter), oq is the Hernquist scale length that 
defines the concentration of matter in NGC 205, and r is 
the distance from the center of NGC 205. Both the mass 
and the Hernquist scale length are variables probed by 
the genetic algorithm. 

The potential of M31 is best represented by a three- 
component m odel containing a bulge, a disk and a halo. 
iGeehan et al.l (2005) derives a "Simple Analytic Bulge- 
Disk-Halo Model" using a spherical symmetric Hernquist 
profile for the bulge (Eqn[2), an infinitesimally thin, ex- 
ponential disk and a spherically symmetric NFW profile 
for the extended dark halo given by 



<^h{r) = -^TiGS^pcrl (y) 



In 



Th 



Th 



(3) 



where 5c is a dimensionless density parameter, pc — 
277.72/i2M0kpc~^ is today's critical density with Hub- 
ble constant h — 0.71 in units of 100 km/s/Mpc, and r ^ 
is the halo scale radius (jNavarro. Frenk fc White 1996^). 
The exponential disk model used by Gechan et aLl 
(|2005f ) is presently without a known analytic solution. 
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Although the equation and its spatial derivatives can eas- 
ily be solved numerically, the amount of time required 
to integrate a set of particles through a complete orbit 
becomes unreasonable for our purposes. An alternative 
approach is to construct a large table of values for the 
spatial derivatives and perform a bi-cubic interpolation 
in both R and z in order to find the force on a single 
particle at each step in the integration. Unfortunately, 
this latter method is unable to significantly decrease the 
amount of computation time required for a single orbit. 
Thus, we use a Miyamoto & Nagai (1975) disk which 
provides results that are comparable to thos e from the 
exponential disk used bv lGeehan et al.l (|2005f ): 



^d{R,z) 



-GMd 



^i?2 + (i?rf + Vi^T&^)2 



(4) 



where Rd is the disk scale length and b is the vertical 
scale factor. 

We use the "Be st-fit Model" values derived by 
iGeehan et al.l (|2005[ ) to describe the various parame- 
ters of M31, with the sole exception of b, the verti- 
cal scale factor, which was not a reported parameter. 
For b we use the vertical scale height of the dust at a 
value of .1 kpc ([H atano et al.' '1997^ . The values re- 
ported by IGeehan e t al. (2005) include the bulge mass, 
with Mb ~ 3.3 X 10^° Mq, the bulge scale factor, with 
o-Q — rh — 0.61 kpc, the disk central surface density, 
with Eq = 4.6 X 10^ Af0 kpc^^, the disk scale radius, 
with Rd = 5.4 kpc, the halo scale radius, with r/j — 8.18 
kpc, and the total mass enclosed inside 125 kpc, with 
M(< 125 kpc)= 5.6 X 10" Mq. To compute the mass of 
the disk we use, 

(5) 



Md = 27ri?5E 



for the halo mass we use, 

Mh = A/(<i25kpc) - Md 



..4x 10^° Mp 







Mb = 4.4 X 10" Me 



(6) 



and to find the dimcnsionless density parameter 6c we 
use the values for r and M^ at r — 125 kpc, and solve 
the mass profile for a NFW halo. 



Mh = AirGdcPcrl 



In 



Th 



r + Th 



(7) 



This results in a dimcnsionless density parameter value 
of, 

^ ^ m 

■iT:pcrl{\n[{r + rh)/rh]-r/{rh + r)) 

= 24.8 X 10''. (8) 

A summary of these parameters is given in Tabled 

3.3. Arrangement of NGC 205 's Test Particles 

In addition to the parameters discussed in Ji l3.1l and the 
potentials outlined in ii l3.2i the satellite's initial, undis- 
torted particle configuration (both in position and veloc- 
ity space) directly affects the appearance of the final, in- 
tegrated distribution of particles. The particles are used 
as tracers of surface brightness and to map the radial 
velocity profile of the galaxy. The present true bright- 
ness profile of NGC 205 lies in between an exponential 
and ri/-* law (|Ke^ [19871: IChoi et al.|[200l . while the 
galaxy's support is thought to come from a mix of both 



rotational and anisotropic velocities (|Geha et al.ll200a) . 
As mentioned in § 12.21 the current observed ratio of the 
maximum rotational velocity to the average velocity dis- 
persion is Wmax/c = 0.21. Hcuce the current structure 
of particles in the satellite lies somewhere between an 
exponential bulge and exponential disk. 

The morphology of NGC 205 before tidal distortion 
is unknown, while a wide range of internal dynamics 
has been observed for cluster dEs, ranging from non - 
rotating to rotational flattened dEs (jGeha et al.l l2003l. 
Furthermore, it is unclear which features of NGC 205 
are due to intrinsic properties, projection, or tidally in- 
teraction. For this reason, we test three configurations 
for NGC 205: a rotating cold exponential disk fii l3.3.ip . a 
non- rotating hot spheroid fii l3.3.3p , and a rotating warm 
exponential disk (§ I3.3.2|) . These configurations explore 
the two extremes, a rotationally supported satellite with 
t^max/c = oo and an isotropically supported satellite with 
Vmax/'y = 0, as well as an intermediate construction, a 
warm disk with 1.05 < Wmax/f < 4.25. These three mod- 
els are outlined in detail below. 

3.3.1. Rotating Cold Disk 

Our most basic model for NGC 205 places mass-less 
test particles in an infinitely thin, fiat disk. The particles 
are randomly distributed in accordance with an expo- 
nential surface brightness profile, as given in Eqn[T] We 
assume that the central (nucleus free) surface brightness, 
/o, has not evolved and that the undistorted NGC 205 
follows an exponential characterized by a variable scale 
length, ro.i. We set the maximum radius of NGC 205 to 
8ro,,;, a location at which very few particles are expected. 
In addition to ro,i, the particle positions in NGC 205 's 
initial disk are also defined by the disk's clockwise rota- 
tion in the x-y plane, 0, and inclination angle, 9. The 
values of these three parameters (rg.i, 0, 0) are obtained 
from the genetic algorithm. 

Initially, the cold disk is completely supported by ro- 
tation, such that fmax/c — oo. Each particle moves in 
a circular orbit about the satellite's center under the in- 
fiuence of a Hcrnquist profile. This results in a circular 
velocity, in the plane of NGC 205's disk, of 



' GMr 
{r + ao)2 



(9) 



where M is the mass of NGC 205, ag is the Hernquist 
scale length of NGC 205, and r is the particle's distance 
from the center of NGC 205. The direction of rotation in 
the disk is determined by the parameter — A6'. Hence, 
the particle's circular velocities are very simply defined 
by their distance from the satellite's center and by the 
three variable parameters: M205, 0205, and 9 — A9. Any 
dispersion found in the final, integrated system is a direct 
result of tidal interaction with M31. 

3.3.2. Rotating Warm Disk 

A more realistic model for NGC 205 arranges mass-less 
test particles in a disk with some vertical thickness. We 
employ the method prescribed by Hernquist ( 1993 ) using 
cylindrical coordinates and distribute the particles in an 
exponentially decaying density profile given by. 



M 



P 



Airrpo 



exp(- 



^i?/ro)sech^ ( — 

,^0 



(10) 
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TABLE 2 
M31 



Parameter 



Symbol 



Value 



Bulge: 



Mass 


A4 


3.3 X 10^0 Mq 


Scale Radius 


n 


0.61 kpc 


Disk: 






Central Surface Density 


So 


4.6 X 10** Mq kpc-2 


Mass 


Md 


8.4 X m^° Mq 


Scale Radius 


Rd 


5.4 kpc 


Vertical Scale Height 


b 


0.1 kpc 


Halo: 






Mass 


Mh 


4.4 X 10" Mq 


Scale Radius 


1-h 


8.18 kpc 


Density Parameter 


&o 


24.8 X 10" 


Total: 






Mass (< 125 kpc) 


*jf(<125 kpc) 


5.6 X 10" Mq 



Note. — M31 paramete r values used in the n umerical simulations 
adopted from the values in lGeehan et al.1 120051 ). 



where R is the cyhndrical radius, M is the mass of the 
sateUite, tq is the disk scale length, and zq is the vertical 
scale thickness. For our simulations, we set vertical scale 
thickness, zq, equal to the scale radius, ro,i, in order to 
create a "puffier" disk. As with the cold disk model 
outlined in Ji l3.3.1i the maximum radius of NGC 205 is 
constrained to 8ro,i and the disk is oriented using the </> 
and 9 parameters obtained from the genetic algorithm. 

The warm disk is supported by a combination of ro- 
tational an d anisotropic ve locities. Using the method 
outlined bv lHernquistI ()1993f ). self consistent disks are ap- 
proximated using the moments of the collisionless Boltz- 
mann equation. In order to avoid imaginary streaming 
velocities (cr? < 0) near the center, the softened version 
of the radial dispersion equation is used. 



to the radial dispersion. 



4 



Aexp(-/R2q^2af/ro) 



(11) 



where a^ is the disk softening radius and A is the nor - 
malization constant. As prescribed bv lHernquistI (|1993D . 
the softening radius is set to ro/4. The normalization 
factor, A, can then be easily computed by specifying a\ 
at i? = ro. 



A 






eM^rl + 2ayro) (12) 



However, we cannot simply set the radial dispersion 
at one scale length equal to the observed dispersion in 
NGC 205 of 35±5km s'^ (|Geha et al.ll2006h . This would 
cause the streaming velocities to become imaginary for 
much of the parameter space. To avoid this potential 
problem, the maximum allowed radial dispersion that 
imposes cr? > is used, with the caveat that crji cannot 
exceed 29 kms^^. The vertical dispersion perpendicular 
to the plane of the disk is given by, 

al = ttG1:{R)zo (13) 

where S(-R) is the disk surface dens ity given by the local 
Hernquist density (|Hernquistlll990D . 



S(i?) = 2pzo = 



MaQZQ 



nR{R + ao) 



(14) 



and ao is the Hernquist scale length. Using the epicyclic 
approximation, the azimuthal dispersion can be related 



^4172 



(15) 



where the epicyclic and angular frequencies are respec- 
tively given by. 



3 a$ 92$ 

'RdR^d^ 



(16) 



(17) 



with the circular velocity Vc given by Eqn [9] and the po- 
tential energy $ of a Hernquist profile given by Eqn [2] If 
we assume that the velocity ellipsoid stays aligned with 
the coordinate axes, the streaming velocity becomes, 

V = ".^ + 4(1 -1^-2^) (18) 

Once these low order moments are determined, the ra- 
dial and vertical velocities for individual particles can 
be drawn randomly from Gaussians with widths cfr and 
az, while the azimuthal velocity is determined by adding 
the streaming velocity, v^, to the velocity value drawn 
from a Gaussian with width a^. This configuration re- 
sults in observable ratios of the maximum rotational ve- 
locity to the average velocity dispersion dispersions of 
1.05 < Wmax/c < 4.25. Since the satellite still contains 
rotational support, the parameter 6 — A9 is used to spec- 
ify the direction of this rotation. 

3.3.3. Non-Rotating Hot Spheroid 

The third model for NGC 205 configures the satellite as 
a dynamically hot, pressure supported spheroid. Hence, 
the velocities of this model are completely isotropic. This 
construction is essentially the counter to the rotation- 
ally supported, cold disk model outlined in § 13.3.11 The 
density profile follo ws the analytical model for spherical 
galaxies outlined bv lHernquistI ()1990D . 



p{r) 



M ao 



1 



(19) 



27r r (r + oq)^ 

where M is the mass of the satellite, ag is the Hernquist 
scale length of NGC 205, and r is the particle's distance 
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from the center of NGC 205. As with the previous two 
models, the maximmii radius of NGC 205 is constrained 
to Sro^i. 

Using the energy dependent distribution function, 
f{E), for a Hernquist profile, a spheroid supported by 
isotropic v elocities can be constructed using random re - 
alizations (jHernquistl Il990t lAscasibar fc BinrievI I2005D . 
The phase-space density as a function of specific energy 



IS, 



fiE) = 



M/al 



47r3(2GA//ao)3/2 



3sin-l(g) + 97(1^^(1 - 2q^){^q^ 



3) 



(l-q2)5/2 

with the dimensionless variable q defined as. 



' GM' 



E 



(20) 



(21) 



The specific energy E oi a particle at position r and 
velocity v is, 

E=-v'^ + <^{r) (22) 

with $(r), the potential energy for a Hernquist profile, 
given by Eqn[21 We normalize the distribution function 
in a volume element dfdv, centered on (r, v), by divid- 
ing fdrdv by its maximum. This maximum corresponds 
to the point in parameter space containing the largest 
amount of mass (and in our case, mass-less test particles). 
Using spherical symmetry, dfdv hecomes Anr'^drAnv'^dv. 
Hence, the maximum of fdfdv is easily found numeri- 
cally and resides at Tpoak = 0.638ao and v = Wc(''poak), 
where Vc is the circular velocity given by Eqn[Hl 

Using the normalized fdfdv, particles are placed in 
phase space using a von Neumann rejection technique 
(jPress et al.l 119861 ) . Tentative values for r and v are 
randomly chosen within the range < r < Sro^i and 

< V < t^esc — ^— 2$(r), and a tentative fdfdv is 
computed. Next, a random number k between and 1 
is selected (where 1 corresponds to the normalized peak 
of fdfdv). The provisional values for r and v are kept 
only if fc < fdfdv, else the process is reinitiated until 
the condition is satisfied. Once r and v are generated, 
position and velocity angular coordinates are randomly 
chosen from < < 27r and — 1 < cos6' < 1, and the 
particle is placed in the satellite. The procedure is then 
repeated until the satellite is completely populated with 
the desired number of particles. This configuration re- 
sults in a galaxy with a density profile given by Eqn 
1191 and that is completely supported by isotropic veloc- 
ities (wmax/c = 0). Note that, unlike the rotating cold 
disk and warm disk, the symmetry of the non-rotating 
spheroid negates the need for parameters </>, 9, and 9—A9, 
while the density profile's dependence on ao reduces ro/s 
use to setting the maximum radius for placing particles. 

3.4. Restricted N-Body Simulations 

The simulations are carried out using a re- 
stricted TV-body approac h (jPfleiderer fc Siedentopjfl96l] : 
iToomre fc Toomrall972[ ) in a Cartesian coordinate sys- 
tem. The use of a Cartesian coordinate system is 
motivated by the format of the available observations: 



IChoi et all (|2002[ ) ph otometric ol p servat ions are defined 
in the x-y plane and iGeha et all ()2006D velocity obser- 
vations are defined in the z direction. We center our co- 
ordinate system on M31 and define the x-y plane as the 
plane of the sky, with y pointing north, along a line of 
increasing declination (DEC), and x pointing east, along 
a line of increasing right ascension (RA). The positive 
z-axis, which is also centered on M31, is therefore along 
our line-of-sight pointing away from us. By this defini- 
tion, an object moving toward us is defined to have a 
negative z-velocity, while an object moving away from 
us is defined to have a positive z-velocity. 

The center of M31 lies at the center of this coordinate 
system (0,0,0). The spherical symmetry of the bulge 
and the halo require no rotation, while the axisymmetry 
of the disk requires a rotation of « = 77.5° around the 
X-axis followed hy a, 9 = 37.7° rotation about the z-axis. 
Thus, if we define a disk as edge on along the cc-axis at 
angles of i = 90° and 9 = 0° then the orientation of the 
M31 disk in the projected plane of the sky is given by: 



X' \ I cos i cos 9 sin 9 

■' ' — ' -cosisin^ cos0 

sin i 





(23) 



For computational purposes, the inverse of the matrix 
above is used to rotate each particle into the frame of 
M31 (with the disk lying in the x-y plane) in order to 
compute the particle acceleration due to the disk, and 
then rotated back into our coordinate system as defined 
above. There are five major steps involved in the simu- 
lation. For each individual orbit, the values of the un- 
known parameters described in S I3.1I are obtained from 
the genetic algorithm as discussed in S I3.5I Each step in 
a simulated orbit is described below: 

1. First, the center of mass (COM) of NGC 205 is in- 
tegrated back in time approximately half an orbit usin g 
a Bulirsch Stoer integration scheme (|Press et al.lll98a l. 
The presence of M31's disk presents difficulties for per- 
forming straightforward calculations of the orbital pe- 
riod. Hence, we approximate the orbital period as twice 
the circular velocity period of NGC 205. Although this 
approximation is far from exact, we intentionally over- 
compensate for the period to ensure that the satellite 
travels at least half an orbit or descends from a rea- 
sonably large radii. If the integration time is longer 
than 6.8 Gyr the orbit is fixed to 6.7 Gyr, since such 
orbits are much longer than the dynamical time scale 
for NGC 205. Since only the COM's x and y positions 
are known precisely, the final z position and the velocity 
vector (vx,Vy,Vz) must be taken from the genetic algo- 
rithm's initial parameter space. The acceleration is com- 
puted using a three-component potential model for M31 
consisting of a bulge, a disk and a halo as described in 
ii l3.2l We set the required accuracy per time step to one 
part in 10^. 

2. At this earlier time, we assume tidal distortion of 
NGC 205 has not yet taken place. After successful in- 
tegration back in time, the initial, undistorted NGC 205 
galaxy is created and populated with 1,000 mass- less test 
particles in a preselected configuration, as discussed in 
S I3.3I The distribution of test particles is not intended to 
be a self consistent tracer of the mass, rather, it is con- 
structed to be a tracer of the expected surface brightness 
profile of NGC 205 prior to its distortion. 
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Fig. 4. — Schematic overview of the simulated interaction between a cold NGC 205 and M31. In Step 1, the center of mass of NGC 205 
is integrated back in time approximately half an orbit. In Step 2, a cold disk representing the undistorted NGC 205 is created. In Step 3, 
the undistorted disk is integrated forward in time in response to NGC 205's and MSl's potentials. In Step 4, a photometric and kinematic 
comparison is done between the integrated system and the observed NGC 205 system. 



3. Next, the undistorted system of particles is inte- 
grated forward in time to the present day under the in- 
fluence of both galaxies' potentials. As discussed in ii l3.2l 
we model NGC 205 as a spherically symmetric static po- 
tential and M31 as a three-component static potential. 
The mass, M205, and Hernquist scale length, 0205, used 
to construct NGC 205 's potential are unknown parame- 
ters and are obtained from the genetic algorithm. We set 
the required accuracy per particle per time step in our 
integration to one part in lO''. 

4. After the integration is complete, the final, simu- 
lated system is compared to the known, observed system. 
The comparison involves a total of 6-7 tests, depending 
on the model used for NGC 205. Although, individually, 
these tests report a probability measurement, their true 
purpose is to provide information to the GA instructing 
it on how to improve subsequent guesses. An outline of 
each test is provided in § 13.61 

5. Steps 1-4 are repeated until the desired number of 
orbits is reached. We then use a genetic algorithm to im- 
prove our initial conditions and repeat the process until 
we can match our observations. For each orbit, param- 
eters resulting in a goodness-of-fit value higher than the 
previous value are recorded. Steps 1-5 are then repeated 
using a new guesses for the free parameters. 

3.5. The Genetic Algorithm (GA) 

Our 10-dimensional parameter space discussed in S I3.1I 
contains approximately 10^^ possible orbits. At an aver- 
age of 6.25 CPU seconds an orbit for our simplest model, 
the cold disk, an iterative search of the space would take 
on order 10^ CPU Gyr. For this reason, we must effi- 
ciently search the parameter space with an optimization 
algorithm. We ha ve chosen to use the optimization tech- 
nique invented bv iHoUaiidl (|1975f l called a genetic algo- 
rithm (GA) for two pragmatic reasons. The first is that 
GAs have the distinct advantage of not getting stuck in 
local optima, a problem that plagues hill climbing tech- 
niques which exploit gradient information. The second 
benefit is that GAs are reasonably effective at reducing 
the amount of searching needed to converge on a solution 



in a given parameter space. Hence, GAs are especially 
useful in problems where iterative searches are compu- 
tational unreasonable. These two features make genetic 
algorithms ideal for searching parameter spaces contain- 
ing an unknown distribution of local optima. 

The inspiration behind the GA is biological evolution. 
The idea follows that evolution is a process determined 
by chromosomes and natural selection. In the GA appli- 
cation, an individual point in parameter space is equiv- 
alent to a chromosome, and the parameters making up 
that chromosome are analogous to genes. As dictated 
by biology, evolution occurs at the point of reproduction 
and individuals with genes advantageous to their species 
are more likely to reproduce. Mutations of genes and 
recombination of chromosomes (crossover) can result in 
offspring that are very different than their biological par- 
ents, and, if these mutations and recombinations result 
in a higher level of success for the offspring, they become 
more likely to reproduce. With each new generation, a 
new population is created and the old population is re- 
placed. 

There are three general modules that define and drive 
all genetic algorithms. These are the Evaluation, Popu- 
lation and Reproduction Modules. Within any particular 
application, the details of each module can vary. The 
purpose of the Evaluation Module is to determine the 
fitness of a particular individual. Fitness is a measure 
of an individual's feasibility as a solution to the problem 
at hand and is commonly determined by a mathemati- 
cal function. The fitness is the only information fed into 
the genetic algorithm and is used in the Population and 
Reproduction Modules. At any given time, the Popula- 
tion Module hosts the current population of individuals, 
called a generation. Encoded in this module are the ex- 
isting population's chromosomes and the methods used 
for creating new generations (e.g. how many children 
are produced by a set of parents). The method of par- 
ent selection (e.g. roulette-wheel, stochastic tournament, 
deterministic sampling) is often contained in this mod- 
ule as well, but can alternatively be a component in the 
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Reproduction Module. The Population Module begins 
by creating an initial population (the first generation), 
which generally consists of random, single-point guesses 
in parameter space. The success of a given GA run only 
weakly depends on this initial first guess. The succes- 
sive replacements of this population by the Reproduction 
Module results in new generations. The Reproduction 
Module determines the method by which new chromo- 
somes will be created. The amount of mutation and the 
method of crossover (the way the parent's chromosomes 
are combined to produce offspring) are defined in this 
module. The resulting offspring becomes an indiv idual 
of the next generation (|Goldberg|ll989l : iDavidllOQll ). 

We use a micro genetic algorithm (micro-GA) driver 
designed bv ICarrolf (|200 1') which employs a tournament 
selection scheme with a shuffling technique to select indi- 
viduals for reproduction. Carroll's micro-GA is charac- 
terized by a small population (5 individuals), with uni- 
form crossover, elitism and no mutation or niching. The 
algorithm works by first arranging individuals in the cur- 
rent population in some arbitrary order (shuffling). Two 
individuals are then randomly chosen from the popula- 
tion, and the fittest amongst the pair is selected as a 
parent (tournament selection). Tournament selection is 
then repeated to select the second parent. Next, genes 
are picked at random from the first or second parent 
to create an offspring (uniform crossover), which subse- 
quently becomes a new member of the next generation. 
The process of selection and reproduction continues until 
the new generation is completely populated. In Carroll's 
micro-GA configuration, elitism ensures that the best in- 
dividual from the previous generation is replicated into 
the new generation and the absence of niching means 
that only one local optima can exist in a population at 
any given time. 



3.6. Measuring Fitness (The Evaluation Module) 

As previously discussed, the only interface between the 
GA and our iV-body simulation is the Evaluation Mod- 
ule, which determines the fitness of an individual. For 
our problem, we use the photometric and kinematic ob- 
servations discussed in §[2] to measure the fitness of each 
simulated orbit. Since the fitness is a solitary number 
fed into the genetic algorithm, each simulated galaxy 
has its photometric and kinematic goodness-of-fit mea- 
surements weighted and combined into a single value. 
The goal of the fitness value is to tell the GA when and 
how to improve. Hence, it is not advantageous to use 
a straightforward x^ method (especially since the simu- 
lated photometric and kinematic results are dependent 
on the same parameters), it is, in fact, more effective 
to weight and combine individual tests that measure the 
qualities of a given orbit. In our simulation, a total of 6-7 
tests are used. After testing many orders and weighting 
schemes we find the following configuration to be most ef- 
fective: two-dimensional two-sided Kolmogorov-Smirnov 
(2DKS) test with a weight of l%(ii l3.6.T|l . slope measure- 
ment with a weight of 9% (S I3.6.2|) . weighted velocity 
profile x^ probability with a weight of 18% (ii l3.6.3p . and 
surface brightness profile x^ probabilities with weights of 
18% each (S I3.6.4|) . The weights and criteria for each of 
these evaluations is discussed below and summarized in 
Tabled 



3.6.1. Two-Dimensional Two-Sided Kolmogorov-Smirnov 
(2DKS) Test 

Our first evaluation uses the two-dimensional two-sided 
Kolmogorov-Smirnov (2DKS) test (Press et al. 1986) . 
The observed surface brightness profile of NGC 205 sys- 
tem is modeled as a distribution of particles placed 
as tracers of the projected surface brightness profile. 
The number of particles in the 'observed' distribution 
is chosen such that it equals the number of particles in 
the simulated distribution. The 2DKS test probability, 
P(D>obs)! is measured using the statistic from Numerical 
Recipes (|Press et al.lll986[ ). 

Simulations resulting a probability > 0.01 are given a 
fitness value of 0.01 and allowed to proceed to the next 
test. Those with values < 0.01 return a fitness value 
equal to the 2DKS probability value. The 2DKS test 
fitness value is thus, 



Fitness — min[P 



(D>obs) 



,0.1] 



(24) 



A weight of only 1% (or a maximum fitness value of 0.01) 
is given to the 2DKS test because it is most effective at 
recognizing bad fits, rather than quantifying good fits. 
That is, the two disadvantages of this test are that it has 
only a 0.20 significance level and a bias towards the center 
(since this is where most of the particles are located) . 

3.6.2. Interior Velocity Slope (Disk Models Only) 

We ensure that NGC 205 has the correct tidally 
undistorted central semi-major axis velocity profile by 
using the velocity slope interior to ~ 270" for the second 
fitness evaluation. This test is performed only on the 
disk models, since the non-rotating spheroid contains no 
preservable internal rotation. To determine the simu- 
lated semi-major axis velocity profile (and measure its 
slope), the minimum distance between each test particle 
and the 5'-shaped major axis of NGC 205 i s computed. 
The slit-masks used by iGeha et all (|2006l ) to observe 
NGC 205's velocity profile are adjacent to the semi-major 
axis, with dimensions 16' (3.8 kpc) long by 4' (1.0 kpc) 
wide. Based on the Geha et al. (2006) collection method, 
simulated particles whose perpendicular distance to the 
major axis is less than 250" (1 kpc) are selected, and 
their velocities are binned and average with neighboring 
particles to create a simulated velocity profile. The ra- 
dial size of each velocity bin along the semi-major axis is 
roughly 1' (0.24 kpc). 

The allowed s lope (m) range is determined from 
IGeha et al.l (|2006[) velocity observations of NGC 205 (see 
? l2.2p . We set a generous slope range for our simulations 
of —0.80 km s^^ arcsec"^ to 0.0 km s~^ arcsec^^ (see 
Figure [3]), corresponding to a —20 km s~^ kpc~^ to 0.0 
km s~^ kpc~^ gradient, if NGC 205 lies at a distance of 
824 kpc. The significance level is chosen such that the 
slope has the maximum range to vary, without ever going 
positive, and corresponds to 1.83cr. 

If the slope falls within the range specified above, it is 
given a probability of 1 and a weight of 9% (or a value 
of 0.09). This is then added to the 2DKS resuh of 0.01 
(for a total fitness value of 0.10) and the code proceeds 
to the next test. If it does not fall within the range, it 
is compared with the slope of a previous run in order to 
determine if the GA is improving and assigns a fitness 
value. Improvement is defined by how the inner slope is 
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TABLE 3 
Determining Fitness, F 



Test 



Passing 
Criteria 



Test's 
Weight 



Fitness if 
Orbit Fails Test 



Fitness if 
Orbit Passes Test 



2DKS Probability 



> 0.01 



1% 



F = P(D > observed)'^ /i = 0.01 



Slope, m 0.0 > m > -0.80 9% 

(in km s"-"^ arcsec"-*^) 



F = /i +0.09(5|- 



nimid 



-r /2 = /i+0.09 



Weighted Velocity 
X^ Probability 



> 10-10 



18% 



F = h 



1.0 



iog[Q(x2k)] 



/3 = /2+ 0.10 + 0.08 xQ(x2|i.) 



Surface Brightness 
X^ Probabilities: 

Radial 
Angular 

Weighted Radial 
Weighted Angular 



> 10-10 

> 10-10 

> 10-10 

> 10-10 



18% 
18% 
18% 
18% 



F = h + f4. + h + h + h 
-1.0 



k 



\og[Qi(x^\v)] 
-1.0 


log[Q5(x"k)] 
-1.0 


loglQeCx^k)] 
-1.0 



iog[Q7(x"k) 



/4 = 0.10 + 0.08 xQ4(x^k) 
/5= 0.10 + 0.08 xQsCx^l!.) 

/6= 0.10 + 0.08 xQ6(x^k) 
/7 = 0.10 + 0.08 xQ7(x^k) 



Where the 2DKS probabihty P(D > observed) is given in iPress et al.| jlQSGl) Used only in the exponential disk models. For the 

0, if|-p---»- 



non-rotating spheroid model, /2 is set to 0.10^ Where nimid — —0.40 km s arcsec and 8 — 



> 1 
m — mmid d 



1, it 



rn„,id 



Where 



< 1 



Qix 1^) is the X probability function from lPress et all (1198611 

changing with respect to previous orbits and is symbol- 
ized by the parameter 6. The parameter S is set equal to 
1 if the absolute value of the slope minus the observation- 
ally best fit slope of —0.40 km s^^ arcsec-^ is less than 
the previous orbit's value (i.e. |mnew — (—0.40 km s-^ 
arcsec"^) I < |mprev — (—0.40 km s~^ arcsec~^)|), and is 
otherwise set equal to 0. The fitness value at this step is 
thus given by, 

0.10, if-0.8< ni < 0.0 kms^i 



Fitness 



0.01-1-0.095 



-0.40 



-0.40) 



otherwise 

(25) 



This results in a minimum fitness value of 0.01 for a 
slope of cxD, and a maximum fitness value of 0.10 for a 
slope very close to 0.0 km s~^ arcsec"^ or 0.80 km s"^ 
arcsec"^. If, instead, the absolute value of the slope mi- 
nus the best fit slope is greater than the previous or- 
bit's value (i.e. |mnew — (—0.40 km s~^ arcsec"^) | > 
liTiprev — (—0.40 km s~^ arcsec~^)| , then the parameter 
6 = and the GA returns a fitness value of 0.01. The 
ultimate goal of this evaluation is to inform the GA if it 
is improving on the internal velocity profile. In the case 
of the non-rotating spheroid, the Fitness is automatically 
set to 0.10 before proceeding to the next test. 

3.6.3. Weighted Semi-Major Axis Radial Velocity Profile 

If the 2DKS test and central slope requirements are 
satisfied, a weighted x^ statistic is determined for the 
semi-major axis velocity profile and a probability Q 



inferred. Similar to the method outlined in ii l3.6.2i 
simulated particles with a perpendicular distance to 
NGC 205's semi-major axis less than 250" (1 kpc) are 
collected and binned. The only difference is that the ra- 
dial length of the velocity bins are slightly larger than 
those used to measure the slope. E ach bin is chosen t o 
be in accordance with those used bv lGeha et al.l (|2006D . 
resulting in an average bin width of approximately 2' 
arcs(M8 kpc). 

The X weights are constructed so that more weight is 
given to the tidally distorted regions beyond ^ 270". 
Ergo, this externally weighted velocity profile compli- 
ments the interior velocity slope discussed in ? I3.6.2I to 
provide a complete description of the semi-major axis 
velocity profile. We use the weighting function, 



w{r) 



■ exp 



ri/2 



:exp 



A 
ri/2 



A 



if r < ri/2 
if r > ri/2 



(26) 



where A determines steepness of the function, and ri/2 
corresponds to the semi-major axis distance at which the 
weight is equal to i. In our model, we set A = 50 and 
^1/2 = 500" in order to produce x^ weights that turn 
on at r^ 270" and rise rapidly. Sin ce the weights are 
discretely sampled in accordance with lGeha et al.l ()2006f ) 
bins, we normalize the weights to the number of available 
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bins, 



A^ X w{ri) 



N 



(27) 



^win) 



where r^ is the semi-major axis value corresponding to 
the center of each velocity bin, and N is the total number 
of bins. These weights are then used to determine the x^ 
statistic, 

^,^^MS^_0^ (28) 

i=i » 

where 5"^ is the average radial velocity found by the simu- 
lation in the ith bin, and both Oi and ct^ are the observed 
average radial velocity and error, respectively, in the ith 
bin. 

Once the x^ value is obtained, its statistical signifi- 
cance can be d etermined using t he x^ probability func- 
tion, Q(x^l^) (iPress et al.lll986f ). where the degrees of 
freedom, v, is given by A'^ — 1. If the weighted radial 
velocity profile probability is > 10"^" it is considered an 
adequate fit and a value of 0. 1-1-0. OSx probability is added 
to the slope and 2DKS fitness'. The code then proceeds 
to the next set of tests. However, if the weighted radial 
velocity profile probability is < 10~^°, then a new fitness 
is computed by subtracting the inverse logarithm of the 
probability from the slope and 2DKS fitness'. So, based 
on the probability, the orbit has an updated fitness value 
of, 



Fitness 



fs = 0.2 



0.10 



^ 0.08 X Qix'^liy) 

1.0 

log[Q(x2|^)]' 



if Q > 10- 



if Q < 10- 



(29) 

This results in a minimum fitness value of 0.10 for a 
Qix^W) probability of 0, and a maximum fitness value 
of 0.28 for a Q{x^W) probability of 1. The contribution 
of the weighted radial velocity profile test to the the total 
fitness is 18%. 

3.6.4. Weighted & Unweigtited Radial & Angular 
Brigtitness Profile 

If the previous three tests (2DKS, slope, weighted ra- 
dial velocity profile) are satisfied, the weighted and un- 
weighted x^ statistics for the projected surface bright- 
ness profile are obtained and probabilities, Q, inferred. 
Since NGC 205 is undergoing tidal distortion, informa- 
tion from both the radial and angular profiles can be 
used as constraints. It is advantageous to look at the 
weighted and unweighted, radial and angular profiles in- 
dependently since they each provide unique information 
about the satellite's orbit. That is, the weighted x^ 
statistics give information about the exterior, tidally dis- 
torted regions of the satellite, whereas the unweighted x^ 
statistics inform about the interior, undistorted regions 
of NGC 205. Likewise, the radial profile reveals the ex- 
tent that the satellite is being compressed or expanded, 
while the angular profile describes the twisting and elon- 
gation of NGC 205. 

In our simulation, mass-less test particles are used 
as tracers of surface brightness. The weighted and un- 
weighted angular and radial distribution of test particles 



is compared with the surface brigh tness profile given in 
Figured which uses the data from lChoi et al.l (|2002f ) in 
the range 100" to 672" and our fitted exponential surface 
brightness profile interior to 100" and exterior to 672". 
In addition to this surface brightness profile, we use the 
positio n angles and ellipticities measured by IChoi et al.l 
(|2002f) from 100" to 672". The position angle and ellip- 
ticities of the isophotes internal to 100" and external to 
672" are held fixed, since interior to this we have removed 
the nucleus and exterior to this the surface brightness 
and sky become comparable. 

For the x^ tests, we step in semi-major axis units of 94" 
(~ 0.38 kpc) and in angular units of 9°. We measure out 
to a maximum semi-major axis distance of 31.2' (~ 7.5 
kpc) and use the center of each areal bin to determine the 
corresponding surface brightness. Since our semi-major 
axis s teps are not coincident with those of IChoi et aLl 
(|2002f ) isophotes, we interpolate to obtain an interme- 
diate ellipticity and position angle that is complemen- 
tary to our step. The unweighted radial surface bright- 
ness profile of the simulated distribution is then found 
by summing the number of tracer particles within each 
predetermined isophotal annulus, while the unweighted 
angular surface brightness profile is measured by total- 
ing the number of tracer particles contained within each 
angular beam. The unweighted x^ statistic for both the 
radial and angular surface brightness profiles can then be 
computed using Eqn[251 setting each weight wt equal to 
1, normalizing our observed distribution Oi to total the 
number of tracer particles used in the simulation, and 
setting af equal to Oi, the number of expected particles 
in the bin. For the weighted x^ tests, the same formula 
is followed, except that the weighting function given by 
Eqn [26] is used to determine Wi . 

Each of the 4 surface brightness profile tests dis- 
cussed above are evaluated independently and combined 
to update the GA fitness. Similar to the velocity pro- 
file test in ii l3.6.31 tests with probabilities > 10~^° are 
considered reasonable fits and are assigned a value of 
0.1 -1-0. 08 X probability. Profiles resulting in probabilities 
< 10^^" are designated bad fits and assigned a value of 
— 1.0/log[(3j(x^|i^)]- The updated fitness then becomes, 

7 

Fitness = /g + ^ /, 

J=4 



0.10 + 0.08 X Qj{x^\v), if Qj > 10 



-10 



-1.0 



iog[Q,{x^w)r 



if Qj < 10-1" 



(30) 



where the j index specifies which of the 4 surface bright- 
ness profile tests, and f^ refers to the weighted velocity 
profile fitness from Eqn [251 Hence, each x^ test can con- 
tribute a maximum of 18% to the total fitness, resulting 
in a maximum fitness value of 1. Once the fitness of a or- 
bit has been determined, the GA can use the information 
to create new, improved orbits. 

4. RESULTS 

In this section, we present the genetic algorithm's re- 
sults for the three NGC 205 models outlined in ^ \'S.'S[ the 
rotating cold disk (§[13]), the rotating warm disk fii l4.2p 
and the non- rotating hot spheroid (ii l4.3p . We perform 
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1000 GA runs on each model (with different initial ran- 
dom number seeds for each run) in order to thoroughly 
explore the parameter space and expose possible degen- 
eracies. Each GA run consists of 1000 generations of 5 in- 
dividual orbits and produces a best orbit (highest fitness) 
from these 5000 orbits. Hence, for each model, 5 x 10^ 
orbits are explored and 10'^ best orbits are returned. Be- 
low, the resulting direction of NGC 205 's approach on 
the sky plane, parameter values and single best orbit are 
given for these three models. 

4.1. Rotating Cold Disk 

We initialize NGC 205 as a rotating, cold exponen- 
tial disk of mass-less test particles, serving as tracers of 
surface brightness. The resulting direction of approach 
(§ 14.1.11) . preferred parameter values (^ I4.1.2p and best 
orbit fi; l4.1.3p for the cold disk simulations are outlined 
below. 

4.1.1. Direction of NGC 205's Approach (Cold Disk) 

Figure \5\ shows histograms of NGC 205's varying di- 
rections of approach on the plane of the sky and selected 
spatial projections for the cold, rotating disk model. The 
histogram in Figure [5^ depicts the best orbit from each 
of the 1000 GA runs. Even without any further con- 
straints, a clear directional preference emerges from the 
1000 orbits. The histogram peaks at 148 orbits with 
NGC 205 advancing from the west (W), followed closely 
by 140 orbits from the west-northwest (WNW) and 130 
orbits from the northwest (NW). Of these 1000 orbits, 
752 are bound to M31, 500 are prograde, 280 are retro- 
grade and 220 are radial (with radial orbits defined by 
arccos|f • v\ < 5°) . 

Figure [Sb shows orbits with additional photometric, 
kinematic, and fitness constraints. The outlined his- 
togram imposes that the weighted velocity x^ proba- 
bility > 10"^*^ and that each of the 4 surface bright- 
ness x^ probabihties return values > 0.1 (ii l3.6.3l -? l3.6.4p . 
This constraint significantly reduces the initial 1000 or- 
bits to 32 orbits and returns a direction of approach that 
lies somewhere between the north-northwest (NNW) and 
west (W), with the peak lying in the NW. Furthermore, 
none of these orbits are retrograde (in fact, none of the 
retrograde orbits have a fitness > 0.25). The enclosed 
shaded region further imposes the condition that fitness 
> 0.86, a value selected to be just below that of the top 
10 orbits (or 1%). This results in a reduction from 32 
orbits to 9, meaning that one of the top 10 orbits fails to 
satisfy all the imposed photometric and kinematic con- 
straints. The remaining directions of approach continue 
to peak in the NW, with only W approaches now ruled 
out. Hence, the directional preference initially suggested 
by the histogram of 1000 orbits is reinforced by the addi- 
tion of photometric, kinematic and fitness requirements. 

More importantly, these constraints rule out other 
possible orbits for t he c old disk model. This includes 
iMcConnachie et aLl (|2004 ) favored 1° long stellar arc-like 
feature (observed in the northwest quadrant of M31) as 
a tidal stellar stream originating from NGC 205. Model- 
ing this trail as tidal debris from the satellite suggests 
encroachment from somewhere in between the north- 
northeast (NNE) and east (E) region (22.5° to 90°) on 
the plane of the sky. While 21% of the orbits from the 
1000 GA runs fall within this 4-bin region, we find that 



these orbits all have poor angular surface brightness dis- 
tributions (i.e they cannot simultaneously match both 
a weighted and unweighted angular x^ probability con- 
straint of 0.1). This result indicates that orbits approach- 
ing from the NNE through E are ruled out as solutions 
by the observed isophotal twisting in NGC 205 (ii l2.ip . 

To further test this conclusion, we perform a single 
cold disk GA run with 3000 particles. While numerous 
GA runs with > 1000 particles are too expensive for our 
purposes, a single GA run allows us to check our results 
against a more complex system. The 3000 particle GA 
run returns a direction of approach from the NNW, an or- 
bit in qualitative agreement with our 1000 particle runs. 
This result is denoted in Figure [5Jd with a dashed line. 

The plots at the bottom of Figure [5] show the trajec- 
tories of the constrained nine orbits corresponding to 
the shaded region in Figure [5)3. All nine orbits have 
very large tangential motions and are not bound to 
M31. Included in the x — y sky projection plot (Figure 
Sfc) are tracers of t h e nin e orbits, an approximation of 



McConnachie et all ()2004D stehar arc, M3rs bulge and 



disk at two scale radii (subtending 5' and 45', respec- 
tively), and two surface brightness isophotes, one at 270" 
corresponding to the onset of NGC 205's tidal distortion 
and the other at 19. '77 indicating NGC 205's tidal ra- 
dius. As discussed, none of the constrained orbits trace 
the stellar arc observed in the northwest quadrant of 
M31. Included in the x — z line-of-sight plot (Figure 
[Hi) are tracers of the nine orbits, M31's bulge and disk 
at two scale radii, and the range of allowed z values (2- 
76 kpc). All nine orbits converge at distances very close 
to M31 and follow roughly similar trajectories. Further- 
more, the nine orbits begin on nearly radial courses with 
0° < /i < 1.5°, where /i is the angle between r and v, the 
position vector from the center of NGC 205 to M31 and 
the velocity vector of NGC 205, respectively. Ergo, the 
radial paths of these nine simulated satellites precludes 
prograde or retrograde encounters. 

4.1.2. Preferred Parameter Values (Cold Disk) 

Figure [6] shows the evolved parameter values resulting 
from the 1000 cold disk GA runs. Each best orbit's pa- 
rameter and corresponding fitness are plotted. Since the 
fitness is determined by combining numerous tests, the 
line illustrating passing criteria for each fitness step is 
also included, with the sole exception of the 2DKS test 
for which all 1000 orbits quantitatively pass. The nine 
constrained orbits from Figure [IJj and the cutoff fitness 
value line for the top 10 orbits are highlighted in red. 
Additionally, the resulting best orbit's parameter values 
from the 3000 particle GA run are included (dashed line) . 

The GA's preference for certain parameter values 
emerges in Figure [6l The most notable concentration is 
in the initial scale radius ro,i of NGC 205, shown in detail 
in Figure [7] (top row). The 1000 simulations converge at 
ro,j = 129 ±12" (0.51 ±0.05 kpc), with orbits passing aU 
fivex^ tests (Fitness > 0.6) at ro,,; = 128±10" (0.51±0.04 
kpc) and the constrained nine orbits (highlighted in red) 
focused at 148 ± 7" (0.59 ± 0.03 kpc), a value equivalent 
to the satellite's present scale radius. This implies that a 
large portion of NGC 205 is unaffected by tides and the 
internal regions of the satellite have experienced little, if 
any, distortion. The remaining parameters experience a 
greater amount of scatter. However, definite trends still 
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Fig. 5. — Cold Disk Model. Direction of approach for NGC 205 modeled as an initially rotating, cold disk, (a) The best orbit from each 
of the 1000 GA runs, (b) Orbits with additional photometric, kinematic and fitness constraints. Included is the resulting best orbit from a 
single GA run with 3000 particles (dashed line), (c) The sky projected [x — y) orbits from the shaded region of the histogram. Also plotted 
are M31's bulge and disk at two scale radii (subtending 5' and 45', respectively), an approximation of McConnachic et al. (2004) stellar 
arc and two surface brightness isophotes, one at 270" (the onset of NGC 205's tidal distortion) and the other at 19. '77 (NGC 205's tidal 
radius). Note, none of these orbits trace the stellar arc observed in MSl's northwest quadrant, (d) The line-of-sight {x — z) orbits from the 
shaded region of the histogram. Also included are MSl's bulge and disk at two scale radii and the investigated z values (2 — 76 kpc). 
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Fig. 6. — Cold Disk Model. Plotted are the parameter values vs. fitness from the 1000 GA runs when NGC 205 is modeled as a rotating, 
cold disk. Each point corresponds to the best orbit from an individual GA run and a particular initial configuration (structure, dynamics, 
orientation) of the progenitor. The included lines denote passing criteria for the various fitness tests: slope test (Fitness > 0.1), weighted 
velocity x^ test (Fitness > 0.2), all five x^ tests (Fitness > 0.6), top 10 (1%) orbits (red line. Fitness > 0.86), and nine constrained orbits 
(highlighted in red). Orbits that fail arc poor fits to photometry and/or kinematics and reside below the Fitness = 0.6 line. Also included 
are the parameter values of the resulting best orbit from a single GA run with 3000 particles (dashed line). 
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Fig. 7. — Cold Disk Model. Shown here are three parameters, for the rotating, cold disk model, that are well constrained by the GA. 
From top to bottom: initial scale length ro,i, velocity magnitude |d| = {v^ + v^ + v^)^'^, and line-of-sight distance z. From left to right: 
parameter values vs. Fitness for the 1000 GA runs, histograms and fitted Gaussians (except z) for all 1000 GA runs, and histograms and 
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kpc. Also included are the parameter values of the resulting best orbit from a single GA run with 3000 particles (dashed line). 



NGC 205 



19 



exist in the data. For example, the simulations favor or- 
bits with very large velocities of 414 ± 94 km s~^ (mid- 
dle row of Figure[71 Fitness > 0.6), preferentially moving 
NGC 205 towards the south {vy = -346 ± 385km s^^) 
and to the east {vx = 333 ± 89km s~^). In addition, the 
GA prefers z distances very close to M31 with z = 11 ± 9 
kpc (bottom row of Figure [71 Fitness > 0.6). An or- 
bit's fitness beyond this distance declines rapidly. Also 
included in Figure [7] are the resulting best orbit's rp.i, 
\v\ and z values from the 3000 particle GA run (dashed 
line). Note that these values are in qualitative agreement 
with the 1000 particle GA runs, with a slightly higher \v\ 
value that is comparable to the velocities of the nine con- 
strained orbits. 

Conversely, some of the parameters cannot be further 
reduced by the simulations. The huge scatter in the ra- 
dial velocity, Vz, indicates that that the GA is unable to 
constrain this value beyond current observations (Figure 
[6l middle row, first column). This result is not surpris- 
ing since the observed error on Vz is quite small. Addi- 
tionally, the GA is unable to reach convergence for the 
satellite's mass, M205, and Hernquist scale length, 0205. 
One explanation for this is that significant degeneracies 
might exist for il/205 and 0205 in the parameter space. A 
second possibility is that perhaps the photometric and 
kinematic fitness tests are not as sensitive to M205 and 
0205 as they are to the other parameters. Hence, the GA 
both produces inconclusive results for M205 and 0205 and 
does not further reduce the observed error on Vz- How- 
ever, it is able to tightly constrain NGC 205's initial scale 
radius and place more general bounds on the remaining 
5 parameters. 

4.1.3. Best Orbit (Cold Disk) 

Figure [5] shows the resulting best orbit from the 1000 
GA runs using the cold disk model. This best orbit 
approaches from the NW, is unbound to M31 (with 



= 546 km s ^ and v^ 



480 km s" 



and has a 



fitness value of 0.90. Figure [5^ displays the simulated 
galaxy and its trajectory on the plane of the sky (in 
black). Also included (in red) are M31's bulge and disk, 
at 2rb and 2Rd, and the pr edicted par t icle d istribution 
for NGC 205, derived from IChoi etall (|2002[ ) isophotes 
and our surface brightness profile (see §[21 Figure[2]). Fig- 
ure [Sb shows this same orbit in the x— z plane, with the 
addition of a vertical red line denoting the allowed z pa- 
rameter values (2 — 76 kpc). Shown in the bottom plots 
are the simulated galaxy's velocity dispersi on (circles) 
and p erformance on the five y^ t ests, with IChoi et all 
(|2002f) isophotes and iGeha et all (|2006[ ) velocity profile 
-|- errors depicted in red. The simulated surface bright- 
ness x^ tests result in probabilities of 0.979, 0.856, 0.999 
and 0.852 for the radial, angular, weighted radial and 
weighted angular tests, respectively, and in a x^ prob- 
ability of 6.0 X 10"^" for the weighted velocity profile. 
Note that while the weighted velocity profile probability 
seems quite low, it satisfies the goal of the simulation 
by matching the turnover and reversal in the semi-major 
axis velocity profile. Furthermore, although this satellite 
began as a purely rotating disk with no dispersion, after 
the simulated interaction with M31 a Gaussian fit to the 
combined velocity distribution for all particles along the 
semi-major axis is 3 km s^^. Even though this simu- 
lated average dispersion is much lower than the observed 



average dispersion of 42 km s^^ (jGeha et al.ll2006D . its 
presence indicates that some of the observed dispersion 
could have been tidally induced by M31. In order to re- 
duce the Poisson noise in the dispersion profile, we reran 
the orbit with 5500 particles. The resulting 5500 par- 
ticle Gaussian fitted velocity dispersion profile is given 
in Figure ^ (solid black line). The parameter values 
of the 1000 particle orbit are M205 = 1.38 x 10^ M^ 



369 km s 



-399 km s" 



52km s-i, 

z = 6.3 kpc, ro,j ="0.6 kpc, </> = 90°, 9 = -66°, and 
^205 — 9.6 kpc. This orbit is projected to pass within 
w 9 kpc of M31's center in the plane of M31's disk. How- 
ever, the radial path of the orbit precludes any sense of 
a prograde or retrograde encounter. 

Using these orbital parameters, we roughly predict the 
position of the satellite and its particles at given points 
in time. Figure [5] illustrates the past, present and future 
predictions for NGC 205 when it is modeled as a cold, ro- 
tating disk. The distribution of particles and positions in 
space are given for times —338 Myr, —169 Myr, present, 
169 Myr and 338 Myr. The panel to the far right in Fig- 
ure[9]shows NGC 205 in the plane of the sky at 338 Myr. 
The simulated satellite experiences a significant amount 
tidal distortion after passing through the disk of M31. 

4.2. Rotating Warm Disk 

This section is identical to §[lll]but for NGC 205 mod- 
eled as a rotating, warm exponential disk of mass-less 
test particles supported by a combination or rotational 
and anisotropic velocities. 

4.2.1. Direction of NGC 205's Approach (Warm Disk) 

In Figure [TOk a distinct directional preference emerges. 
The histogram peaks at 228 orbits with NGC 205 advanc- 
ing from the northwest (NW), followed by 151 orbits from 
the north- northwest (NNW) and 89 orbits from the west- 
northwest ( WNW) . Almost half the orbits are contained 
within these three directional bins. Furthermore, these 
directions of approach reinforce those from the 1000 cold 
disk model runs, which peaked in the west (W). Of these 
1000 orbits, 808 are bound to M31, 233 are prograde, 593 
are retrograde and 174 are radial. 

Figure [TOb shows an effective reduction of 1000 or- 
bits to 5 and returns a direction of approach lying some- 
where between the north (N) and northwest (NW) , with 
the peaks in both the NW and NNW. Of these 5 orbits, 
2 are prograde, 1 is retrograde, and 2 are radial. The 
shaded region enclosed within the histogram further im- 
poses that the Fitness > 0.76, a value selected to be just 
below that of the top 10 orbits (or 1%). This constraint 
eliminates the retrograde orbit and reduces the 5 orbits 
to 3, indicating that 7 of the top 10 orbits fail to satisfy 
all the imposed photometric and kinematic constraints 
given above. The remaining directions of approach con- 
tinue to peak in the NNW, with only N approaches now 
ruled out, a result similar to that of the cold disk model's. 
Hence, the initial directional preference from the 1000 
orbits histogram is reinforced with the addition of pho- 
tometric, kinematic and fitness constraints. 

These restrictions rule out other possible orbits, includ- 
ing those tracing the stellar arc-like feature seen to the 
north of M31. The warm disk model simulations place 
15% of the 1000 GA orbits in bins NNE through E. How- 
ever, only 2 of these 153 orbits meet both a weighted and 
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Fig. 10. — Warm Disk Model. Same as Figure[5]but for NGC 205 modeled as a warm disk with both rotational and anisotropic velocities. 
Note that the directions of approach on the plane of the sky are similar to those from the cold disk models. 



22 



Howley et al. 



unweighted angular x^ constraint of 0.1 and fail to meet 
a weighted and unweighted radial x^ constraint of 0.1. 
As with the cold disk model, orbits tracing the observed 
stellar arc are ruled out as solutions by NGC 205 's pho- 
tometric profile. 

4.2.2. Preferred Parameter Values (Warm Disk) 

Preferences by the GA for various parameter values 
can be seen in Figure [11] As with the cold disk model, 
the most notable focus is in the initial scale radius VQi 
of NGC 205, shown in detail in Figure [T^] (top row). 
Both the 1000 simulations and the orbits passing all five 
X^ tests (fitness > 0.6) converge at rg,,; — 108 ± 11" 
(0.43±0.04 kpc), with the constrained three orbits (high- 
hghted in red) focused at 113 ± 13" (0.45 ± 0.05 kpc). 
These GA selected scale radii are close to the satellite's 
present scale radius of 148" (0.59 kpc), which imphes 
that the internal regions of the satellite have experienced 
very little distortion. There is a greater amount of scat- 
ter in the remaining parameters, however, distinct trends 
still exist in the data. Shown in the middle row of Figure 
[T2l orbits with very large velocities of 368 ± 78 km s~^ 
(Fitness > 0.6) are favored, moving NGC 205 primar- 
ily towards the southeast {v^ = 229 ± 54 km s^"'^ and 
Vy — — 312± 132 km s^^). In addition, the GA prefers z 
distances very close to M31 with z = 11 ± 8 kpc (bottom 
row of Figure [121 Fitness > 0.6). Beyond this distance 
an orbit's fitness declines rapidly. However, some pa- 
rameters are not well constrained by the simulations. As 
with the cold disk model, the radial velocity v^ contains 
a significant amount of scatter amongst the allowed pa- 
rameter values. In addition, the satellite's mass, M205, 
and Hernquist scale length, 0205, are also unresolved pa- 
rameters. Hence, as with the cold disk results, the GA 
both produces inconclusive results for M205 and 0205 and 
does not further reduce the observed error on Vz- How- 
ever, it is able to tightly constrain NGC 205 's initial scale 
radius and place more general bounds on the remaining 
five parameters, which also happen to closely match the 
values from the cold disk simulations. 

4.2.3. Best Orbit (Warm Disk) 

Figure [13] displays the resulting best orbit from the 
1000 GA runs for the warm disk model. This or- 
bit approaches M31 from the NNW, is unbound (with 
\v\ = 514km s~^ and Wcsc = 475 km s~^), and has a 
fitness of 0.87. The simulated surface brightness x^ tests 
result in probabilities of 0.730, 0.994, 0.999 and 0.657 
for the radial, angular, weighted radial and weighted 
angular tests, respectively, and in a x^ probability of 
5.0 X 10~^° for the weighted velocity profile. This satel- 
lite is initialized with a radial dispersion at one scale 
length of (Tfj-ra i = 3.7km s^^ and a vertical dispersion 
of crz=ro i = 5.5 km s~^. After the simulated interaction 
with M31, the semi-major axis velocity profile contains 
an average observable dispersion of 6 km s~^ (which is 
slightly greater than its initial dispersion), with a maxi- 
mum dispersion of 4 km s^^ at a radius of 2.1 kpc. The 
parameter values of this orbit are A/205 = 1.0 x 10^ Mq, 



269 km s" 



-435 km s 



52 km s 



z = 7.2 kpc, ro,,; = 0.5 kpc, = 83°, 6* = -77°, and 
1205 — 7.5 kpc. As with the cold disk model, this orbit 
is projected to pass within w 9 kpc of M31's center in 



the plane of M31's disk. Also, the radial path of this 
orbit precludes any sense of a prograde or retrograde en- 
counter. 

Figure [T4l illustrates the past, present and future pre- 
dictions for NGC 205 when it is modeled as a warm disk 
supported by rotation and anisotropic velocities. The 
distribution of particles and positions in space are given 
for times -338 Myr, -169 Myr, present, 169 Myr and 
338 Myr. The panel to the far right in Figure [Ml shows 
NGC 205 in the plane of the sky at 338 Myr. The sim- 
ulated satellite experiences a significant amount of tidal 
distortion after passing through the disk of M31 . 

4.3. Non-Rotating Hot Spheroid 

This section is similar to Ji l4.1l fc S I4. 21 except that here 
NGC 205 is modeled as a non-rotating, hot spheroid of 
mass-less test particles supported by isotropic velocities. 
Since this configuration of particles has increased veloc- 
ity dispersion (compared to the disk models), we include 
and additional discussion about our ability to match this 
observed component in § 14.3.41 

4.3.1. Direction of NGC 205's Approach (Hot Spheroid) 

In Figure [T5I a clear directional preference is seen in the 
initial histogram of 1000 orbits. The histogram peaks 
at 432 orbits with NGC 205 advancing from the north- 
northwest (NNW), followed by 178 orbits from the north 
(N) . More than half the orbits are contained within these 
two directional bins. Furthermore, these directions of ap- 
proach reinforce those from the disk model runs, which 
peaked in the W (cold) and NW (warm). Of these 1000 
orbits, 653 are bound to M31 and 308 are on radial tra- 
jectories. 

The top-right panel of Figure [TS] shows a histogram 
of orbits with added photometric, kinematic, and fit- 
ness constraints. The outlined histogram imposes that 
weighted surface brightness x^ probabilities return val- 
ues > 0.1, and that the unweighted surface brightnes 
and weighted velocity x^ probabilities > 10^^°. No- 
tice that these constraints differ from those applied to 
the exponential disk models. Since the hot spheroidal 
model follows a Hernquist profile characterized by an r^' ** 
law (and not an exponential profile), we relax the sur- 
face brightness x^ conditions that are more sensitivite 
to the satellite's inner regions, a portion that has ex- 
perienced little, if any, tidal distortion. These imposed 
constraints reduce 1000 orbits to 5 and return a direction 
of approach lying somewhere between the north (N) and 
north- northwest (NNW), with a peak of 4 orbits in the 
NNW. The shaded region enclosed within the histogram 
further imposes that the Fitness > 0.685, a value selected 
to be just below that of the top 10 orbits (or 1%). This 
constraint reduces the 5 orbits to 3, indicating that 7 of 
the top 10 orbits fail to satisfy all the imposed photomet- 
ric and kinematic constraints given above. The remain- 
ing directions of approach all fall in the NNW, ruling out 
N approaches, a result similar to the disk models' find- 
ings. Hence, the initial directional preference suggested 
by the 1000 orbits histogram is reinforced with the addi- 
tion of photometric, kinematic and fitness constraints. 

These restrictions rule out other possible orbits, includ- 
ing those tracing the stellar arc-like feature seen to the 
north of M31. The hot spheroid simulations place a mere 
2.5% of the 1000 GA orbits in the 4 bins NNE through 
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Fig. 11. — Warm Disk Model. Same as Figure[6]but for the warm disk model with only three constrained orbits, highlighted in red. 



E. None of these orbits meet a weighted angular x^ con- 
straint of 0.1. Hence, orbits tracing the observed steUar 
arc are ruled out as solutions by the twisting of elliptical 
isophotes in the tidally distorted regions of NGC 205. 

4.3.2. Preferred Parameter Values (Hot Spheroid) 

The resultant parameter values from the 1000 hot 
spheroid GA runs are shown in Figure 1161 The spher- 
ical nature and density profile of this model effectively 
reduces the number of parameters from 9 to 6, eliminat- 



ing ro,i, (/), and 9. 

The GA's preference for certain parameter values is 
quantified in Figure 1171 The top two rows (Fitness 
> 0.6) show the well resolved tangential components of 
the velocity vector with Vx = 163 ± 77 km s~^ and 
vy = -376 ± 46 km s-\ effectively moving NGC 205 
towards the southeast. This results in a net velocity of 
\v\ = 417 ± 33 km s'^ (Fitness > 0.6). Contrastly, only 
the magnitude of the velocity could be well resolved in 
the case of the cold disk model. The GA also prefers z 
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Fig. 12. — Warm Disk Model. Same as Figure [71 but for the warm disk model. Orbits with Fitness > 0.6 give ro,i = 109 ± 11", 
\v\ = 368 it 79 km s^-"^, 2 = 11 ± 8 kpc, values nearly identical to those from the cold disk models. 
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Fig. 13. — Warm Disk Model. Same as Figure[8]but for the warm disk model whose resulting best orbit has a fitness of 0.87. 
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Fig. 14. — Warm Disk Model. Same as Figure [9] but for the warm disk model. 
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Fig. 15. — Hot Spheroid Model. Similar to Figures [Sl and [Tol but for NGC 205 modeled as a non-rotating hot spheroid. For the spheroid 
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model the unweighted photometric 

the directions of approach on the plane of the sky are similar to those found for the disk models. 
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distances very close to M31 with z = 9 ± 5 kpc (bot- 
tom row of Figure \T7\ Fitness > 0.6). As with the disk 
models, an orbit's fitness beyond this distance declines 
rapidly. In the bottom-right plot of Figure I16[ a corre- 
lation of higher fitness to lower 0205 values can be seen. 
This is a result of the initial particle distribution's depen- 
dence on 0205 in Eqn[Tni a dependence that did not exist 
for the disk models. Hence, the GA's attempt to repro- 
duce the exponential surface brightness profile results in 
this correlation. However, Figure \T^ also demonstrates 
that some of the parameters are not well constrained by 
the simulations. As with the disk models, the radial ve- 
locity Vz contains a significant amount of scatter and is 
likely not constrainable beyond observations. In addi- 
tion, the satellite's mass, M205 is unresolved. Hence, as 
with the disk results, the GA produces inconclusive re- 
sults for M205 and does not further reduce the observed 
error on Vz- However, it is able to reasonable constrain 
NGC 205's tangential velocity and line-of-sight distance, 
values which also happen to closely match those from the 
disk simulations. 

4.3.3. Best Orbit (Hot Spheroid) 

Figure [TH] displays the hot spheroid's resulting best 
orbit from the 1000 GA runs. This orbit approaches 
M31 from the NNW, is just barely bound (with |i;| = 
434km s~^ and Uosc = 483 km s~^), and has a Fitness 
of 0.73. Notice that on the bottom-right plot the inner 
velocity profile is hatched out indicating that the slope 
test was not used for this model. The simulated surface 
brightness x^ tests result in probabilities of 8 x 10^^°, 
6x 10"^", 0.757 and 0.83 for the radial, angular, weighted 
radial and weighted angular tests, respectively, and in 
a x^ probability of 0.002 for the weighted velocity pro- 
file. The orbit's poor performance on the unweighted 
surface brightness x^ tests is a direct result of compar- 
ing the simulated spheroid, constructed using a Hern- 
quist density profile, to the observed NGC 205, which 
follows an exponential profile. Despite this fact, the or- 
bit performs quite well on the weighted x^ tests that 
model the tidally distorted regions of NGC 205. After 
the simulated interaction with M31, the semi-major axis 
velocity profile contains an average observable dispersion 
of 17 km s~^ (which is lower than the observed disper- 
sion of 42 km s~^), with a maximum dispersion of 36 
km s~^ at a radius of 1.0 kpc. The parameter values of 
this orbit arc M2Q5 = 1.5 x 10^ M©, v^ = 168km s"\ 



= -397km s- 



= 57 km s ^, 



5.8 kpc, and 



1205 — 1-7 kpc. As with both disk models, this orbit is 
projected to pass within « 9 kpc of M31's center in the 
plane of M31's disk. Since this satellite is not rotating, 
prograde or retrograde definitions do not apply. 

Figure fT9l illustrates the past, present and future pre- 
dictions for NGC 205 when it is modeled as a hot 
spheroid supported by isotropic velocities. The distribu- 
tion of particles are given for times —384 Myr, present, 
76 Myr and 384 Myr. The gravitationally bound satel- 
lite traces out a loop orbit. The panel to the far right 
in Figure [12] shows NGC 205 in the plane of the sky at 
384 Myr. The simulated satellite experiences a signifi- 
cant amount of tidal distortion after passing through the 
disk of M31 . 

4.3.4. Semi-Major Axis Velocity Dispersion 



Figures [H [13] and [18] show the best fit orbit for each 
of the three models. However, none these orbits pro- 
duce profiles that match the observed velocity dispersion 
profile. This is expected since our fitness tests do not 
contain a velocity dispersion profile test, meaning our 
orbits are not encouraged to match this observable. The 
main reason this feature is not used is because the low 
number of particles (1, 000) used for each orbit makes it 
nearly impossible to fit a histogram and recover an accu- 
rate measurement for the velocity dispersion at all points 
along the semi-major axis of NGC 205. 

We expect that of the three models the hot spheroid 
model, supported completely by isotropic velocities, is 
most capable of reproducing the observed dispersion. To 
test this, we use 3, 000 mass-less test particles plus the 
best fit hot spheroid model's velocity and line of sight 
distance parameters given in ? I4.3.31 and we conduct a 
coarse, iterative search in the mass (M205) and Hern- 
quist scale length (0205) dimensions. We search these 
two parameters because they are unconstrained in all 
three models and heavily infiuence the simulated veloc- 
ity dispersion profile. After producing discrete relations 
between at and the mean velocity dispersion, we fit a 
power law line and interpolate to find the values of ah, 
at each value of M205, that result in a mean dispersion 
of 35 km s~^. This result is displayed as circles in the 
top panel of Figure [20] The solid line in the top panel 
of Figure [50] indicates the analytic solution for these pa- 
rameter values in a tidally undistorted satellite with a 
projected velocity dispersion of 35 km s~^ at a radius of 
1 kpc (Hcrnquist (1990), Eqn [42]). Although this fine 
measures the velocity dispersion at a given radius, rather 
than the mean velocity dispersion, the result suggests 
that tidal interactions between at hot spheroid progen- 
itor and parent galaxy have not significantly impacted 
the velocity dispersion within a 1 kpc radius. 

In addition to matching the observed mean dispersion, 
we find that all these combinations of M205 and ah suffi- 
ciently reproduce the velocity dispersion profile, includ- 
ing the observed dip at the center. We randomly select 
one of these points (A/205 — 4.2 x 10^ Mq, ah — 1.2 kpc), 
denote it with a black circle and display the orbit in the 
bottom 4 panels of Figure [20] The poor fitness values 
reported by this orbit are not concerning since this orbit 
was not a true result of the optimization scheme. Hence, 
while the mass and Hernquist scale length parameters are 
not constrained in the optimized parameter space search, 
their relative values are suggested by the observed veloc- 
ity dispersion. 

5. DISCUSSION 

This paper presents the methods and results of an op- 
timized, restricted A^-body search for NGC 205 's orbit 
and system parameters. We account for uncertainties in 
the parameter space by placing very liberal upper and 
lower limits on our parameters, using three dynamically 
distinct models for the initial configuration of NGC 205, 
and running our optimization algorithm (GA) 1000 times 
on each model. We find that out of the 10^^ possible or- 
bits, not only are certain orbits and portions of parame- 
ter space carved out, but this convergence is also model 
independent. These findings are outlined below: 

• The simulations indicate that NGC 205 is ap- 
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Fig. 18. — Hot Spheroid Model. Same as Figures [Sl and [131 but for the hot spheroid model whose resulting best orbit has a fitness of 0.73. 
(d) Note that the inner velocity profile is hatched out indicating that it was not used to determine the orbit's fitness. 
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Fig. 19. — Hot Spheroid Model. Same as Figures [51 and [TH but for the hot spheroid model. The positions are given at times —384 Myr, 
present, 76 Myr and 384 Myr. In this model, NGC 205 traces out a loop orbit. The panel on the right is a zoomed in view of NGC 205 at 
384 Myr on the plane of the sky. 
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Fig. 20. — Hot Spheroid Model. The circles in the top panel show, for the tidally distorted best-fit hot spheroid orbit, the combinations 
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preaching from the NW region on the plane of the 
sky, moving towards the SE (i.e. increasing in RA, 
decreasing in DEC), with the cold disk model fa- 
voring an approach from the NW and the warm 
disk & hot spheroid models favoring the NNW. 

• These orbits do not trace out nor do they intersect 
the 1° long stellar arc observed to north of M31, 
which was hypothesize d as a tidal stream emanat - 
ing from NGC 205 bv iMcConnachie et all ([200l . 
Orbits that trace out this region on the plane of the 
sky (i.e. approaches from the NNE through E re- 
gion) poorly match the photometry of NGC 205, 
specifically, the isophotal twisting. Instead, ap- 
proaches from the NW-NNW region are preferred, 
thus ruling out the stellar arc as a trail of debris 
from NGC 205. 

• Large velocities in the range 300-500 km s~^ are 
highly favored for NGC 205. While the major- 
ity of the orbits found by the GA are bound, 
the best fitting orbits are either very close to es- 
cape velocity or unbound. However, such veloci- 
ties are not completely uncommon for local group 
satellites. Proper motion measurements reveal 
that the Large Magellanic Cloud is very near es- 
cape velocity a nd likely on its first passage about 
the M ilky Way (jKallivavalil et a l. 2006; Bcsl a et all 
l2007f ). The dwarf spheroidal galaxies And XIV 
and And XII are near, or exceed, their local es- 
cape speeds and are presuma bly falling into the Lo- 
cal Group for the firs t time (jMaiewski et al.l[2007t 
IChapman et al.ll2007t l. The large tangential com- 
ponent of NGC 205 's simulated velocities translates 
into a proper motion of ~ 0.1 mas yr~^. 

• The GA is able to better fit NGC 205 's observed 
kinematic profile when the satellite's line-of-sight 
distance is very close to M31. The resulting 
fits report relative distances between NGC 205 
and M31 in the range 2-20 kpc. Note that 2 
kpc is our lower bound on the line sight dis- 
tance parameter. These dista nces correspond to a 
change of O.OTTtaois mag in IMcConnachie et all 
()2005l ) measurements of NGC 205 's red giant 
branch tip (TRGB), analogous to either an in- 
crease in their adopted TRGB absolute magni- 
tude (M™^^) and/or reddening correction {Aj), 
or in a decrease in the TRGB observed magni- 
tude (/trgb)- Furthermore, this magnitude dif- 
ference is in a greement with the j = 0.07 mag er- 
rors quoted bv IMcConnachie et al.l (|2005l ) on the 
individual distance measurements to NGC 205 and 
M31. As an aside, a magnitude difference could 
also be accounte d for by NGC 205's comp osite stel- 
lar population (jSalaris fc Girardil [2005) resulting 
from a recent (^ 0.5 Gyr) burst of star formation 
(jMarleau et al.ll2006h . Given these possible correc- 
tions and the close relative distances found by the 
simulation, the question is raised: Could NGC 205 
he in front of M31 rather than behind it? 

• The disk models return tightly constrained initial 
scale lengths in the range 0.39-0.55 kpc, values 
approximately equal to NGC 205 's current scale 



length of 0.59 kpc. This result is consistent with 
the photometric and kinematic observations which 
suggest the central region of NGC 205 is largely un- 
affected by tides. Hence, both the inner kinematics 
and photometry are determined by the initial pa- 
rameters of the system. 

• The cold disk model demonstrates its plausibility 
as a configuration for NGC 205 by the induction of 
velocity dispersions into its initially cold profile via 
tidal interactions. Hence, all three models contain 
some dispersion along their final semi-major axis 
velocity profile. 

• For the best fit hot spheroid model, we can re- 
produce the observed mean velocity dispersion of 
35 km s~"'^and the semi-major axis velocity disper- 
sion profile by using specific combinations of two 
otherwise unconstrained parameters, satellite mass 
(Af205) E^nd Hernquist scale length (a205). 

• For all three models, NGC 205's large tangential 
velocities and close approaches to M31 indicate 
that it is difficult to distort a satellite to the ex- 
tent that NGC 205 is observed. In order to repro- 
duce the abrupt turnover and reversal in the veloc- 
ity profile, NGC 205 must come very near to M31. 
While the photometric profile can be reproduced 
at large distances, the resulting velocity turnover is 
less abrupt and the reversal is not as steep. How- 
ever, at very close approaches (< 20 kpc), it is very 
easy to disrupt NGC 205. In order to avoid com- 
plete disruption and/or long tidal streamers, the 
satellite must approach M31 with a very high ve- 
locity. Hence, the combined photometric and kine- 
matic observations have managed to remove degen- 
eracies in some of the parameter space. 

• The constrained orbits are primarily radial and 
projected to pass within 10 kpc of M31's center. 
For the cold disk model, orbits that are not radial 
and have a fitness > 0.25 are prograde. In contrast, 
a larger fraction of the warm disk model's orbits 
with fitness > 0.25 are retrograde. This result is 
not surprising since the warm disk model contains 
some initial velocity dispersion, allowing its satel- 
lites to be more easily distorted. Thus, some of the 
warm disk satellites are able to reproduced milder 
photometric dist ortions by ap proaching on retro- 
grade orbits (Re ad et al.ll2006l ). 

• Given that many of the best fitting orbits are 
both radial and unbound, it is not surprising that 
NGC 205 has only experienced a small amount 
of distortion. That is, an unbound orbit indi- 
cates that NGC 205 is on its first passage while 
satellites on radial orbits are more diffi cult to dis- 
rupt than those on prograde orbits (|Read et al.l 
l2006f ). Although the future projections show 
NGC 205 whizzing through M31, inclusion of dy- 
namical friction would change this result by slow- 
ing t he satellite after hit ting M31 for the first time 
(Seguin fc Du Draall994 ). This braking mechanism 
is so efficient that it is likely NGC 205 will eventu- 
ally merge with M31. 
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Although the GA's search was quite successful, it was 
unable to constrain the mass and dark matter content. 
This latter result could be due to either degeneracies in 
the parameter space or an insensitivity to these param- 
eters by the fitness tests. Furthermore, it appears that 
to some degree the algorithm did get stuck in local op- 
tima, a result evidenced by the scatter in the solutions 
and the fact that some orbits had very low fitness. Al- 
though we attempt to circumvent this undesirable result 
by running the GA numerous times, we cannot be 100% 
certain that we converge on the "best orbit" in the pa- 
rameter space. However, we do find that other GA and 
fitness test configurations result in declining levels of a 
best orbit fitness' and are thus confident that our method 
is effective. While other optimization tools are available 
(i.e. simulated annealing), they too have their problems 
and limitations, and cannot offer guaranteed improve- 
ments. 

6. FUTURE DIRECTIONS 

We have presented the results of optimized restricted 
A''-body simulations of NGC 205 's tidal interaction with 
M31 using initially self-consistent satellite models and 
parameter spaces tailored specifically to the system. Al- 
though this simple model produces satisfactory results, 
improvements to the simulations will generate a more 
detailed and accurate account of the interaction, possi- 
bly produce tighter constraints on the parameters, and 
provide greater insight into the past, present and future 
of NGC 205. In this final section we discuss two ob- 
servational advances that can be made and a handful 
of numerical improvements that can be applied to the 
NGC 205-M31 system in order to better and more accu- 
rately model their interaction. 

The simulations can be significantly improved by re- 
ducing the observational errors. As discussed, the 
current value for the line-of-sight distance between 
NGC 205 and M31 is e stimated to be 39 ± 37 kpc 
(|McConnachie et al.ll2005l ). An improvement on the rel- 
ative distance between NGC 205 and M31 would sig- 
nificantly constrain the parameter space. In addition, 
if NGC 205 indeed has a large tangential motion of 
300 — 500 km s~^, the satellite's proper motion could be 
resolved at long wavelengths over the course of a decade 
or two, thus providing another significant constraint on 
the system. These two measurements, in conjunction 
with current sky positions, radial velocities, and mod- 
els of M31's mass distribution, would provide a complete 
description of NGC 205 's orbit. 

Our simple model of NGC 205 's interaction with M31 



yields results that are consistent and relatively indepen- 
dent of model type. An obvious improvement to the 
current simulation is to use velocity dispersion as a con- 
straint on the system. Another improvement, applying to 
only the spheroid model, would be to weight the bright- 
ness of the mass-less test particles such that the initial 
surface brightness profile matches the observed exponen- 
tial profile of NGC 205. Furthermore, this system can 
be better modeled with the inclusion of dynamical fric- 
tion and mass loss, effects which act to slow the system 
down and increase the level of distortion in the satellite. 
The addition of these processes effectively transforms our 
time-independent model into a time-dependent problem. 
Another modification is to construct NGC 205 as a self 
gravitating model with both dark and luminous parti- 
cles, thus switching from a restricted A^-body to a full 
A^-body simulation. However, this more realistic model 
changes the problem from 0{N) to a 0{N'^), leading to 
a notable increase in computation time and making it ex- 
tremely costly to run the simulation along with a genetic 
algorithm. Despite this complication, a self-gravitating 
model has a significant advantage in that it allows for 
evolution in NGC 205. This construction provides the 
necessary platform to test various formation scenarios 
for the satellite since it can be initialized as a spiral or 
dirr and then transformed into a dE through tidal inter- 
actions. If this improvement were to be made, one might 
also consider adding an initial gas mass to the progeni- 
tor galaxy and performing hydrodynamic simulations for 
completeness. Last but not least, NGC 205's orbit can be 
improved by adding the perturbations from M31's other 
satellites, such as M32, as well as the potential of the 
Milky Way 
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